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ABSTRACT 

Recent studies have shown that the cross-correlation coefficient between galaxies and 
dark matter is very close to unity on scales outside a few virial radii of galaxy halos, 
independent of the details of how galaxies populate dark matter halos. This find- 
ing makes it possible to determine the dark matter clustering from measurements of 
galaxy-galaxy weak lensing and galaxy clustering. We present new cosmological pa- 
rameter constraints based on large-scale measurements of spectroscopic galaxy sam- 
ples from the Slo an Digital Sky Survey (SDSS) Data Release 7 (DR7). We generalise 
the approach of iBaldauf et al.l (|201Cll ) to remove small scale information (below 2 
and 4/i~^Mpc for lensing and clustering measurements, respectively), where the cross- 
correlation coefficient differs from unity. We derive constraints for three galaxy samples 
covering 7131 deg2, containing 69150, 62150, and 35088 galaxies with mean redshifts 
of 0.11, 0.28, and 0.40. We clearly detect scale-dependent galaxy bias for the more 
luminous galaxy samples, at a level consistent with theoretical expectations. When we 
vary both erg and fim with other cosmological parameters fixed (and marginalise over 
non-linear galaxy bias), the best-constrained quantity is cr8(f^m/0.25)°'^^ = 0.80±0.05 
(1(7, stat. -|- sys.), where statistical and systematic errors (photometric redshift and 
shear calibration) have comparable contributions. These strong constraints on the dark 
matter clustering suggest that this method is competitive with cosmic shear in current 
data, while having very complementary and in some ways less serious systematics. We 
therefore expect that this method will play a prominent role in future weak lensing 
surveys. When we combine these data with WMAP7 CMB data, constraints on erg, 
fim, w^c and '^rrii, become 30-80 per cent tighter than with CMB data alone, 
since our data break several parameter degeneracies. 

Key words: gravitational lensing: weak - cosmology: observations - cosmological 
parameters - large-scale structure of Universe - cosmological parameters 



1 INTRODUCTION 

The currently accepted cosmological model that is broadly 
consistent with multiple observations, known as ACDM, 
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is dominated by dark ingredients: dark matter, which we 
observe through its gravitational eflects, and dark energy, 
the presence of which was inferred due to the accelerated 
expansion of the universe as detected using supernovae 
iRiess et al.lll998l : |Perlmutter et al.lll999l '). Further attempts 
to constrain this model, such as those described by the 
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Dark Energy Task Force jAlbrecht et al.ll2006l '). rely on ob- 
servational methods that can broadly be classified in two 
ways: geometric measurements such as supernovae (standard 
candles) and Baryon Acoustic Oscillations (BAO, standard 
rulers); and measurements of large-scale structure growth. 
The latter measurements of structure growth - particularly 
as a function of time - can constrain the initial amplitude of 
matter fluctuations, the matter density, and even the nature 
of dark energy; the scale-dependence of structure growth can 
be used to constrain the neutrino mass. 

Theoretical predictions for structure growth, such as 
from perturbation theory or A''-body simulations, are clean- 
est when expressed in terms of fluctuations in the den- 
sity of dark matter. Fortunately, weak gravitational lens- 
ing provides us with a way of observing the total mat- 
ter density (including dark matter), via the deflections 
of light due to intervening matter along the line-of- 
sight, which bot h magnifles and distorts galaxy shapes 
(for a review, see Bartelmann fc SchneideJ 200 ll : lRefregie"rl 
I2OO3I : iHoekstra fc Jainll2008l : iMassev et al.ll201(]| ). The lens- 
ing measurement that is commonly used to constrain 
the amplitude and growth of matter fluctuations is 'cos- 
mic shear', the auto-correlation of galaxy shape distor- 
tions due to intervening matter along the line-of-sight. 
Since the initial detections of cosmic shear a decade ago 
llBacon et al.ll2000l:IVan Waerbeke et al.ll200(]| : [Rhodes et al] 
I2OOII : iHoekstra et al.ll2002l '). increasingly enlarging datasets 
and sophisticated measurement techniques have led to 
stead ily decreasing errors, both statistical and systematic 
(e.g., IFu et aklliooa : ISchrabback et al.ll201Gl ). 

However, cosmic shear is, by its very nature, a difficult 
measurement: in the auto-correlation of galaxy shape distor- 
tions, coherent systematic errors (such as those induced by 
seeing or distortions in the telescope) become an additional 
additive term. Moreover, intrinsic alignments with the local 
densit y field that anti-corre late with the real gravitational 
shear (|Hirata fc Seliakll2004l ) can contaminate cosmic shear 
measurements in ways that are difficult to remove. 

iBaldauf et al.l ||201Q ) provided an alternate approach to 
constraining the growth of structure using gravitational lens- 
ing which is less subject to the aforementioned difficulties. 
This approach involves the combination of two measure- 
ments: the auto-correlation of galaxy positions (galaxy clus- 
tering), and the cross-correlation between foreground galaxy 
positions and background galaxy shears (galaxy-galaxy lens- 
ing, which measures the galaxy-mass cross-correlation). By 
combining these two measurements, we can recover the mat- 
ter correlation function, the quantity that is most easily 
predicted by the theory. To reduce uncertainties associ- 
ated with exactly ho w galaxies populate dark matter halos, 
iBaldauf et al.l (|201G| ) construct a two-point observable that 
explicitly eliminates all information below scales equal to a 
few times the typical dark matter halo virial radius. The use 
of these two observations allows for a direct measurement 
of the galaxy bias (the factor relating the matter and the 
galaxy density fiuctuations, which can be both mass- and 
scale-dependent), thus eliminating one of the main system- 
atic uncertainties in using galaxy clustering alone to con- 
strain the matter power spectrum, by converting it to a 
statistical error over which we marginalise when constrain- 
ing cosmology. This measurement can constrain the ampli- 
tude of matter fluctuations at quite low redshift, which is 



very useful when combined with higher-redshift measure- 
ments, providing a measure of structure growth in the time 
when dark energy is most dominant. Also, since it relies on 
shear cross-correlations rather than auto-correlations, coher- 
ent additive errors in galaxy shapes can be removed from the 
analysis entirely. 

This paper i s a p roof of concept of the method described 
in lBaldauf et al.l (|2()1Q ) to constrain the amplitude of matter 
fluctuations at 2 < 0.4, using data from the Sloan Digital 
Sky Survey (SDSS). For this measurement, we use lens sam- 
ples that have spectroscopy, primarily Luminous Red Galax- 
ies (LRGs) due to their homogeneous photometric proper- 
ties, simple selection criteria, and the large cosmological vol- 
ume that they sample. By dividing our sample into three lens 
samples, we can test for consistency between the results at 
different redshifts (modulo the expected amount of evolution 
due to the different mean redshifts). We will demonstrate 
that even this very shallow survey can constrain the ampli- 
tude of matter fluctuations at the ~ 6 per cent level, which 
is especially cosmologically interesting when combined with 
Cosmic Microwave Background (CMB) data. 

We begin in Sec. [2] with a more detailed outline of the 
theoretical background behind the observation we wish to 
carry out, and simulations that we use for tests of this 
method. The data that we use is described in Sec. [S] and 
our observational methodology in Sec. |4] The observational 
results for the galaxy-galaxy lensing are in Sec. [S] and for 
the galaxy clustering, in Sec. [G] We show the resulting con- 
straints on cosmological parameters and on galaxy bias in 
Sec.[71 and conclude in Sec. [H] with perspective on how this 
method may be used in upcoming surveys that will carry 
out deep, wide-field lensing observations. 

Here we note the cosmological model and units used in 
this paper. All estimates of observed quantities assume a 
fiat ACDM universe with n„i = 0.25, = 0.75; we dis- 
cuss the implications of this choice in Sec. 12.3.31 Distances 
quoted for transverse lens-source separation are comoving 
/i~^Mpc, where Ho = 100 /i kms"^ Mpc"^ Likewise, AE 
is computed using the expression for in comoving co- 
ordinates, Eq. ([7)1. In the units used, Hq scales out of ev- 
erything, so our results are independent of this quantity. 
Finally, 2-dimensional separations are indicated with capi- 
tal R, 3-dimensional radii with lower-case r (occasionally r 
may denote r-band magnitude as well; this should be clear 
from context). 



2 THEORY 

The most basic theory predictions for the growth of struc- 
ture are phrased in terms of the statistics of the matter dis- 
tribution - for example, the 2-point matter auto-correlation 
function S^-m-mir) or the power spectrum Pmm(fc). Here the 
matter auto-correlation function is defined in terms of the 
matter density contrast 5m = pm/pm — 1 as 



Cmm(r) = (5m(a;)5m(a; -|- r)). 



(1) 



Perturbation theory is sufficient to predict such statistics of 
the matter distribution when the perturbations are linear 
(density contrast (5m ^ 1); A^-body simula tions are used to 
predi ct the non- linear power spectrum fe.g.. lHeitmann et al.l 
l2010i) in the absence of modifications due to gas physics, 
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which may be significant on the scales use d for typi- 
cal weak lensing analyses (Zhan & Knox 2004; Jing ct al.' 
20061 : iRudd et alBoOS : . Zentnc r ct al, 2008: Semboloni et al.. 
201ll ). 

Galaxy redshift surveys allow us to constrain analogous 
auto-correlation functions for the galaxy density field, Cgg(^) 
or Pgg(fc). Unfortunately, the connection between the theory 
predictions for the matter statistics to the two-point statis- 
tics of the galaxy density field is non-trivial. We can define 
the relation as 



^gsir) = b (r)Cnim(r). 



(2) 



On large scales, it is possible to use the linear bias approx- 
imation, 6(r) — constant, where the galaxy bias depends 
on the mass of the dark matter halos hosting the galaxies. 
However, the bias is also scale-dependent on smaller scale s 
jCole et al.ll2"005l : ISmith et al.ll2007l : ISanchez fc Coi3l2008ll . 
< 20/i~^Mpc for galax;ies in very massive halos. The exis- 
tence of galaxy bias causes significant difficulty in inferring 
the statistics of the underlying matter density field from 
galaxy redshift surveys, additional information is needed. 

Galaxy-gala:xy weak lensing provides a simple way to 
probe the connection between galaxies and matter via their 
cross-correlation function 

Cgm(r) = (5g(a;)C(a; + r)). (3) 

This cross-correlation can be related to the projectecfl sur- 
face density around lensing galaxies 



i?2 + n' 



dn, 



(4) 



where 11 is the line-of-sight separation measured from the 
lens, and therefore = R^ + Il^. This surface density is then 
related to the observable quantity for lensing, the tangential 
shear distortion 74 of the shapes of background galaxies, via 



where 



AE(i?) = 7t(i?)E, = S(< 7?) - E(i?), 



— 2 r^' 

E(< i?) = — / i?'di?'E(i?') 

R Jo 



(5) 



(6) 



When averaging over ('stacking') large numbers of lens 
galaxies to determine the average signal around them, the 
resulting matter distribution is axisymmetric about the line- 
of-sight. The observable quantity AE can be expressed as 
the product of two factors, a tangential shear 74 and a geo- 
metric factor 



Ec 



Da{Zs 



(7) 



AnG Da{zi)Da{zi, Zs){1 + zlY 

where Da{zi) and Da^Zs) are angular diameter distances 
to the lens and source, Da{zi,Zs) is the angular diameter 
distance between the lens and source, and the factor of (1 + 
zl)~^ arises due to our use of comoving coordinates. 

Generally, for some 2-point statistic ( (for example, 
the real-space correlation function ^ or Fourier-space power 
spectrum P(fc)), we can relate the three possible 2-point 



^ In Eq. ^ we ignore tlie radial lensing window, which is so broad 
as to be insignificant on all but tfie largest scales, as was demon- 
strate d explicitly in the context of this method by iBaldauf et al.l 
ll2010l 'l. 



correlations that can be constructed out of the matter and 
galaxy fields, ^mm, Cggi ^^i^d Cgm, as follows: 



6(C) 



Cgg — Cmm — Cgm- 



(8) 
(9) 



All quantities in these equations are a function of scale, 
where the scale depends on the exact statistic (e.g., 3D r, 
2D R, Fourier wavenumber k, multipole £). Here is the 
galaxy bias relating the galaxy and dark matter fluctua- 
tions, and ric\ defined as Tcc^ — Cgm/"^ CmrnCggi the CrOSS- 

correlation coefficient between the matter and galaxy fluc- 
tuationfl Generically, the galaxy bias tends to a constant 
value on large scales ('linear bias'), and the cross-correlation 
coefficient approaches one, but the rate at which this hap- 
pens depends on the choice of statistic ^. In particular, if 
^ is defined as a product of either a Fourier mode (i.e. the 
power spectrum) or of a count in cell (of varying size, called 
the smoothing size), then Ircc'l < 1 (note that no shot-noise 
subtraction is applied here). In this case, the deviation of 
from unity can be related to stochasticity, which is de- 



„(C) 



fined as 



{{S^-b^^^5^f) = Cgg-26(«'Cgm + (&'^')\„ 



(10) 



This is zero if ri'p — 1. However, the rate at which rii^ ap- 
proaches unity as a function of either the size of the cell or 
the wavevector of the Fourier mode is slow, because stochas- 
ticity (such as the shot noise caused by finite number of 
galaxies) contributes to it. This rate of convergence to unity 
is even worse if compensated windows with positive and neg- 
ative weights, such as for the aperture mass statistic, are 
used (|Schneider et al.ll 19981): this eff ect h as been observed i n 
practice bv. e.g.. lSimon et all (|2009l ') and I Julio et~aLl l|2012l ). 

On the other hand, ^ can be defined as a correlation 
function, as in Eq. (O, not as a product of a field with it- 
self (or another field) , in which case the shot noise does not 
explicitly contribute to it except at zero lag. A related statis- 
tic in Fourier space is the shot-noise-subtracted power spec- 
trum, where stoc hasticity is ex p licitly subtracted. In this 
case, as shown in iBaldauf et al.l (|201G| ). ri'p is much closer 
to unity (except at zero lag) and the scale dependence of 
is significantly reduced (which is why shot-noise subtraction 
is a standard procedure in the analysis of the galaxy power 
spectrum). Moreover, even on scales where fe'''^ is strongly 
scale dependent, ri'p is close to unity, with deviations from 
unity of only a few per cent on scales above 3/i~^Mpc, where 
the scale-dependent bias can be tens of per cent. In this case, 
rcP has no relation to stochasticity, since its contribution 
does not enter or is explicitly subtracted from it, and we no 
longer need to have Ircc'l < 1. 

If we can ensure that we are working in a regime where 
the cross-correlation rcP « 1, or, more generally, if we have a 
robust model for its scale dependence, then we can infer the 
combination of the mean matter density and the correlation 
statistic of matter by combining the galaxy-galaxy lensing 



^ This statistic is often denoted r. We use the subscript 'cc' to 
avoid confusion with 3D length scales. 
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and galaxy clustering statistics, 

^=Pm[r<^f Cm^. (11) 

So far, t his discussion has been fairly general. 
iBaldauf et all (|2010l ) carried out a detailed exploration of 
the behaviour of rcc' for a variety of statistics ^, using a 
simulated sample of Luminous Red Galaxies residing in dark 
matter halos with masses > 3 x 10^^ Mq, at z = 0.23. As 
shown there, a key point in determining the optimal statistic 
is that we want to avoid information from within the halo 
virial radius, because those are the scales for which the corre- 
lation coefficient is intrinsically quite different from unity in 
a way that cannot be predicted from first principles (without 
some detailed model for how galaxies populate dark-matter 
halos). The observed lensing signal AE is therefore quite 
non-optimal from the perspective of wanting to do cosmol- 
ogy using large scales only, because as shown in Eqs. (O 
and ([6}, at a given R it depends on the surface density of 
matter around galaxies all the way fro m R = 0. 

The statistic that w a s prop osed bv lBaldauf et al.1 (I2OI0I ) 
and iMandelbaum et al.l l|2010l ) to remove small-scale infor- 
mation is known as the annular differential surface density 
(ADSD) T, defined as 

r{R;Ro) = AE(i?)-(^:^^ AE(i?o) (12) 

= ly^'^'^iR') (13) 
[R^EiR) ^ RlE{Ro)] . 

This statistic depends not only on projected separation R, 
but also on some scale _Ro; as demonstrated in Eq. (|13p . 
T{R;Ro) is completely lacking any information from below 
Ro . We thus have to choose a value of Ro that is appropriate 
for our particular application. We will consider several Ro 
values in this paper, but generally we would like this to be a 
few times the typical dark matter halo virial radius (a point 
that we exa mine in more d e tail in Sec. l2.4p . As demonstrated 
in detail in iBaldauf et all (|20ld ). the advantages of such a 
choice are that (a) the correlation coefficient ri'^^ ~ 1 for all 
scales R ^ Ro, and (b) the few per cent deviations from 1 
can be calculated quite accurately via perturbation theory 
(which is only applicable in this regime outside of halo virial 
radii). It was shown that the deviations of rcc from unity 
can be described well with one free parameter related to 
non-linearity of the bias, 62. In this paper, we allow the data 
(specifically galaxy auto-correlations) to determine 62 , which 
will in turn determine the small deviations of rcc from unity. 
In addition, because T is a partially compensated statistic, it 
is not very susceptible to issues that can plague the projected 
correlation function (uigg) such as sampling variance from 
large-scale modes uniformly shifting Wgg up or down. 

The approach described here, which entails removing 
small-scale information completely, is a conservative ap- 
proach that minimises systematic uncertainties due to all 
the things we do not know on small scales (how galaxies 
populate dark matter halos, baryonic effects on the matter 
power spectrum, etc.) at the expense of increasing the sta- 
tistical errors. Alternative approaches involving halo occu- 
pation distribution (HOD) modeling have also been consid- 



ered (|%^t^ ] |2006l : ICacciato et"aLll2009l : ILeauthaud et al.l 
I2OIII , |2012| )" as ways to combine the lensing and cluster- 
ing observations to constrain cosmology. Those approaches 
can potentially give smaller statistical errors, since they use 
the small-scale lensing signals which typically have the best 
S/N, but they are subject to additional systematic uncer- 
tainties both in terms of theory interpretation and obser- 
vational uncertainties that are more pronounced on small 
scales (e.g., intrinsic alignments; magnification; data pro- 
cessing challenges near bright lens galaxies). 

To be more quantitative, our approach is that the small- 
scale galaxy auto-correlation contains no cosmological infor- 
mation, since there is nothing in the distribution of galaxies 
within the halo that has a simple relation to cosmological pa- 
rameters. While small scale galaxy clustering can constrain 
HOD models, this by itself does not help in cosmological 
constraints. Moreover, it is potentially dangerous to rely on 
small scale clustering to constrain cosmological models, be- 
cause one can never be sure that the HOD parametrisation 
is sufRciently general and that there is no artificial breaking 
of degeneracies with cosmological parameters due to insuf- 
ficient generality. HOD models explored to date may not 
allow all possible freedom in terms of how can galaxies be 
placed inside the halos. Moreover, they ignore issues such as 
assembly bias (explicit dependence of clustering properties 
on assembly history, rathe r than just mass alone; iGao et ali 
I2OO5I , iGao fc Wh"it^l2007l ') that can change the relation be- 
tween small- and large-scale clustering information. Once we 
are trying to place cosmological constraints at the ~ 5 per 
cent level, where these small details become more important, 
it is more robust simply to remove the small scale clustering 
regime. Our approach explicitly does that. 

When testing our procedure, we will apply it to a real- 
istic simulated mock sample, which we have generated using 
an HOD model known to reproduce the galaxy two-point 
correlation function, but our claim is that our procedure 
should work on any sample. 

While we believe that the galaxy auto-correlation con- 
tains no useful information on small scales, the galaxy-dark 
matter correlation does contain information on halo mass, 
which in combination with the galaxy auto-correlation on 
large scales can provide i ndependent c o smolo gical informa- 
tion using the method of ISeliak et al.l (|2005 '). Our current 
method cannot take advantage of this additional informa- 
tion from the small scale lensing, so in this sense it is sub- 
optimal. 

One additional aspect to our approach that is meant to 
reduce systematic uncertainties is that we do not simply use 
all of our lens galaxies in one large sample to get a small 
statistical error. Instead, we have several lens samples at 
different redshifts. This way, we can check for consistency of 
the results with each sample, and check for deviations from 
our assumptions about rcc or observational systematics that 
scale with redshift (such as our understanding of the source 
redshift distribution, which is more important when the lens 
redshift approaches the typical source redshifts). 

2.1 Simulations 

While we argued in the previous section that our approach 
is, by design, fairly insensitive to the details of how galaxies 
occupy dark matter halos, it is nevertheless useful to test 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy-galaxy lensing cosmology 5 







h 


o-g rig «)dc 


0.25 


0.75 


0.7 


0.8 1.0 -1 



Table 1. Cosmological parameters adopted for the simulations: 
matter density relative to the critical density, dark energy density 
parameter, dimensionloss Hubble parameter, matter power spec- 
trum normalisation, primordial power spectrum slope, and dark 
energy equation of state p = w^^p. 

the whole procedure on a mock data sample that is as close 
as possible to the real data. Here, we repeat the descrip- 
tion of the A'^-body s imulations that were used for valida- 
tion of the method in iBaldauf et all ()20ld ) and that we use 
for additional tests in this paper. We use the Ziirich horizon 
'zHORIZON' simulations, a suite of forty pure dissipat ionless 
dark matter simulations of the ACDM cosmology (jSmithl 
I2OO9I ) . Each simulation models the dark matter density field 
in a box of length L = 1500/i"^Mpc, using A^p = 750^ dark 
matter particles with a mass of Mp = 5.55 X 10"/i"^Mq. 
The cosmological parameters for the simulations in Table [T] 
are inspired by the resul ts of the WMAP cosmic m icrowave 
background experiment (ISpergel et al.ll2003l . 120071 ). The ini- 
tial co nditions were set u p at redshift z — 50 using the 2LPT 
code (IScoccimarrd [l99^ . The evolution of the A^p equal 
mass particles under gravity was then foll owed using th e 
publicly available TV-body code GADGET-2 (|Springel|[2005h . 
Finally, gravitationally-bound structures were identified in 
each simu lation snapshot using a Friends-of-Friends (FoF, 
iDavis et' al. 1985) algorithm with linking length of 0.2 times 
the mean inter-particle spacing. We rejected halos contain- 
ing fewer than twenty particles, and identified the potential 
minimum of the particle distribution associated with the 
halo as the halo centre. In total, we identify halos in the 
mass range 1.1 x lO"/i"^M0 ^ Mvir ^ 4 x 10^^/i~^Mq. 

We populate the halos in these simulations with galaxies 
using the Halo Occupation Distribution (HOD). This model 
requires us to specify probability distributions for (a) the 
number of galaxies in our sample that occupy a halo of masfl 
M, and (b) the radial distribution of galaxies within halos. 
The HOD can be separated into terms representing 'central' 
galaxies living at the centre of halos, and 'satellite' galaxies 
that are distributed more widely within the halofl We as- 
sume that a halo can only contain a satellite if it also has 
a central galaxy. This assumption may not be entirely valid 
for a colour-selected sample such as LRGs, if the central 
galaxy is very bright but slightly too blue to be included in 
the sample. This will have effects on scales below the virial 
radius: the galaxy-dark matter correlations will be reduced 
on very small scales. It will also reduce galaxy clustering 
in cases when this satellite has another satellite in the same 
halo. We expect these effects to become small on scales larger 
than the virial radius. 

Details of the five-parameter HOD that we used, and 

^ In principle the number of satellite galaxies could depend on 
other parameters such as formation time; however, the HOD does 
not include dependence on anything other than mass. 

Due to limited resolution, we do not attempt to place the 
satellites in subhalos, but rather distribute them probabilistically 
within the host halo. 



tests of how well it describes the samp le abundance, l e nsing , 
and clustering statistic s, are given in iBaldauf" et all (|2010[ ) 
and iReves et al.l (|2010l ). The satellite fractions range from 
10 per cent at the lower luminosity end, to ~ 5 per cent at 
higher luminosity. These results are c onsistent wi th previous 
estimates of LRG environments fe.g.. lReid fc Spergcl 2003 ). 



2.2 Non-linear bias model 

The analysis in IBaldauf et all l|2010l ) was focused on mod- 
eling the cross-correlation coefficient r^c, since this is the 
only quantity that is needed to relate the measurement of 
galaxy auto-correlation and galaxy-dark matter correlation 
to the dark matter auto-correlation. However, it is useful to 
analyse how well we model clustering and g-g lensing data 
separately with a given non-linear bias model. One reason 
to do so is that this allows us to choose different minimum 
scales (-Ro) for galaxy-galaxy lensing and galaxy clustering. 
We expect that lensing information will be quite insensitive 
to the details of HOD modeling: both a satellite and a cen- 
tral galaxy give approximately the same g-g lensing signal 
for separations larger than the virial radius. So, we would ex- 
pect the lensing signal to be fairly model-independent down 
to the virial radius. In contrast, the clustering signal will 
depend very sensitively on how satellite galaxies are popu- 
lated within the virial radius, so the clustering signal up to at 
least twice the virial radius will be quite model-dependent. 
For example, if there are no satellites, then the clustering 
signal drops to almost zero within twice the virial radius, 
while if all the halos have one central galaxy and one satel- 
lite radially distributed as the dark matter, then the cluster- 
ing signal is similar to the lensing signal. This sensitivity to 
how satellite galaxies populate halos suggests that we should 
choose a larger value of Ro for clustering than for lensing 
measurements. A second reason to use larger Ro for clus- 
tering is that the statistical errors are significantly smaller 
than for lensing, so the analysis of clustering data is much 
more sensitive to inaccuracies in the theoretical model. A 
third reason is that if we can model both of these functions 
with a few free parameters, we can use the better-measured 
galaxy clustering data to determine these parameters. 

To do so, we require models not just for rcc, but also 
for scale-dependent bias, that are as realistic as possible. In 
order to interpret the measurements without taking ratios of 
noisy quantities, we must have some well-motivated way of 
describing the non-linear bias of the sam ples that we study. 
We consider the same loca l bias m odel of lFrv fc Gaztanagal 
( .1993. ) as in IBaldauf et all (|2010l ). 5h = fe5m + 62^^ + hSL 
which contains a local bias relation between galaxy and dark 
matter density up to third order and combines it with stan- 
dard perturbation theory (SPT), which expands the den- 
sity perturbation into a series (5m = ^m' + ^m' + where 
^m'' is the Gaussian linear theory prediction and is of 
order [5m']", to calculate the next-to-leading order correc- 
tions to the galaxy-galaxy and galaxy-matter power spectra. 
The third order bias enters only through a renormalisation 
of the leading order bias parameter, and does not have an 
explicit influence on the observable correlators. At the next- 
to-leading order, the corrections to the galaxy-galaxy and 
galaxy-matter power spectra come from the auto-correlation 
of (5m and its cross-correlation with the second order density 
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1.1 r 

1.08^ 




Figure 1. Top: Correlation function cross-correlation coefficient 
between the mock LRGs and the dark matter (data points), with 
errors determined from the box to box variations in simulations. 
Lines are predictions from our model using different values of 62 • 
Middle: Tgm for mock LRGs and the model predictions. Bottom: 
Tgg for mock LRGs and the model predictions. 

(2) 

erturbation (5m . In evaluating these terms, we can define 



pertu 
JSmil 



Smith ct al. 2009) 



A{k) = 



(27r)3 

B{k)^j-^,P{q)P{\k-q 



F2iq,k-q)P{q)P{\k-q\) (14) 

(15) 



where F2 is the SPT mode coupling kernel (see, e.g., 
iBernardeau et al.ll2002l ) 

5 2 
-F2(<7i, Qa) = y "(<7i> <72) + y/3(<7i, <72)> (16) 

with 

(17) 

Upon Fourier transforming, we obtain the correspond- 
ing correlation functions, which are given by 

izz{r,z) = &^Cmm,NL(r,2:) + 2662?.4(r,z) + ib2Cs(7-,2) (18) 



and 



^gm(r-, 2:) = 6Cmm,NL(r, 2:) + 62CA(r, z). 



(19) 



^,a(t) and ^B{r) are the Fourier transforms of A{k) and 
B{k). In principle, ^mm.NL should be the correlation func- 
tion corresponding to one loop perturbation theory. Taking 
SPT at face value, the Fourier transform is ill behaved and 
we replace it by the non-linear correlation function measured 
in the A'^-body simulations. Note that = Ciin- While there 
is a single &2 in Eqs. 1)18^ and (|19p . we find a better fit to the 
simulated data if we have two separate values for the two 
equations, 62, gg and 62, gm that are related via 



^'2,gm = 62 



0.2. 



(20) 



We will explore the implications of this choice on the real 
data, as well. 

As shown in iBaldauf et all (|2010l ). the above model 
can be used to predict the cross-correlation coefficient in 
the linear and weakly non-linear regime. It predicts rcc to 
be unity on large scales and to drop below unity as one 
goes to smaller scales, with explicit functional form given 
by rcc' = 1 — {h2/b)'^(,i\-n{r). We know this model to be im- 
perfect in the sense that other non-linear bias parameters 
at a quadratic and cubic lev el may be needed to properly 
model the data (|Chan et al.i r2012: Baldauf ct al. 2012), but 
these higher order parameters may not be that different in 
terms of their effect on the scale dependence of the statistics 
we study here, so we group them into a single parameter 62 
for the purpose of this paper. 

As seen in the top of Fig. [TJ this model describes the cor- 
relation coefficient down to 3/i~^Mpc, below which physics 
from within the virial radius begins to affect the results. As 
argued above, we expect that these effects are more signifi- 
cant for the auto-correlation than for the cross-correlation. 
We must choose the minimum scale at which we can still 
adopt this model. Our method for doing so will be described 
in Sec. 12.41 it is based on carrying out our analysis on sim- 
ulated data, and checking that we can recover the true cos- 
mology in the simulations. Before we can do so, we next 
describe how we model the observable quantities in real and 
simulated data, T^e and T^m- 



2.3 Modeling the observables 



Our approach is to use Eqs. (|18p and (|19|) to model the 
two observables. The data are used to constrain linear fe, 
quadratic bias 62 , and the dark matter power spectrum times 
the matter density, as in Eq. (|11[) . This is the full, non- linear 
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matter power spectrum, as shown in iBaldauf et all (|2010t ). 
We will use Monte Carlo Markov Chain (MCMC) methods, 
in which the data are compared to the model, hence for each 
set of cosmological parameters we must compute the fully 
non-linear dark matter power spectrum. 

2. 3. 1 Matter power spectrum 

We obtain the estimated linear ^mm(»') by specifying the cos- 
mological paramet ers using the CAMB linear gravity solver 
(jLewis et al.ll2000l ). which is part of the cosmomc package 
that we use for the estima tion of the cosmological parame- 
ters iLewis fc Bridlell2002l ). We increase the accuracy of the 
solver, by setting accuracy_level=l . 5, and checked that in- 
creasing accuracy_level does not change our results. The 
correlation functions ^nim('~) are calculated at the effective 
redshifts of the three galaxy samples, given in Table [21 

To compute the non-linear power spectrum for arbitrary 
cosmological parameters that differ from our fiducial ones 
(00 ), there are several possible approaches that we could 
take (given our simulations that are on a grid of cosmo- 
logical parameters) . The change in cosmological parameters 
affects the non-linear power spectrum in two ways: first, via 
the change in the linear power spectrum; and second, via 
the change in non-linearity corrections. Since the first ef- 
fect is dominant, we account for it accurately using analytic 
calculations of the linear matter power spectrum, only in- 
terpolating on our simulation grid to account for the second 
(much smaller) effect. If we relate the non-linear and linear 
correlation functions via 

U{r\e) = Ciin(r|0)|4lSr = Cun{r\0)a{r\0), (21) 

then we can Taylor expand a{r) around our fiducial cosmo- 
logical parameters. 



Uirm =5iin(r|0) 



a{r\0o) + Y^ 



da{r\0) 



A9, 



- ^un{r\0)a{r\0o) 



E 



d log a{r\0) 



A9, 



(22) 



where the index i runs over the parameters for which we 
have ^mm on a grid {as, ris, 0,m., Hq, and redshift z). As an 
example of how this works for one parameter, changing ug 
mostly changes the amplitude of the correlation functions 
by the square of the ratio of two values of erg under con- 
sideration and this change is propagated exactly. Only the 
second-order change in the shape of the non-linear correc- 
tions is Taylor-expanded. 

Our fiducial model has Qm ~ 0.25, erg = 0.8, ris — 1.0, 
h = 0.7, and z — 0.23; we have 8 simulations of this cos- 
mology. To obtain the derivatives of non-linear corrections 
with respect to cosmological parameters (needed in Eq. l22|) . 
we use further models with Qm ~ (0.2,0.3), erg = (0.7,0.9), 
Us = (0.95,1.05), Ho = (65,75) km s^^ Mpc-\ and 7 dif- 
ferent redshift slices between z — and z — 0.51. Non-linear 
correlation functi ons for these models are obtained from A'^- 
body simulations (|Smithll2009l '). For each non-fiducial model, 
all parameters but one are kept at the fiducial value. For the 
parameters for which we have two simulations bracketing the 
fiducial value {Qm, erg, z, Hq), we use different derivatives 
depending on whether the corresponding value for the target 



model is above or below the fiducial value. We opted to do 
this rather than calculating the second derivative to avoid 
numerical errors blowing up when the quadratic term be- 
comes dominant. By construction, such modelling exactly 
reproduces the non-linear matter correlation function for 
models for which we have simulations. 



2.3.2 Massive neutrinos 

We would also like to place constraints on massive neutri- 
nos, which requires some additional corrections to the for- 
malism in Sec. 12.3.11 We parametrise the neutrino mass ef- 
fect as the sum of masses for the three neutrino families, 
m„, and include three different ways that they affect the 
matter power spectrum. First, lensing is sensitive to the 
total gravitational potential, which includes a contribution 
from massive neutrinos. This requires us to use the Pois- 
son equation to relate potential to density perturbations, 
the latter of which must include the neutrino contribution 
{5p — Pcdm^cdm + Pb^b + Pi/^i/, whcrc cdm, b and u sub- 
scripts denote cold dark matter, baryons and neutrinos, re- 
spectively). For = 0.15eV, f„ = p^/(pcdm + Pb) = 0.6 
per cent, and since neutrino perturbations go from 5^ = Scdm 
on large scales to (5^ = on small scales, this effect sup- 
presses the weak lensing power spectrum on small scales by 
1.2 per cent. 

The second effect is the usual suppression of matter 
fluctuations due to the fact that neutrino fluctuations are 
suppressed on small scales, which in turn leads to a sup- 
pressed growth of cold dark matter and baryon fluctuations. 
For = 0.15eV, this effect leads to a 8 per cent suppres- 

sion of the matter power spectrum. The two effects combined 
thus lead to 9.2 per cent suppression. 

The third effect is the non-linear evolution correction, 
which further enhances this effect. For k < 0.1/i/Mpc the 
effect of neutrinos on the matter power spectrum in the lin- 
ear regime ca n be described as a reduction of the amplitude 
and a red tilt (|Bird et al.ll2012l '). For = 0.15eV, this is 

a 4 per cent reduction in erg and -0.01 reduction in ris. Most 
of the mode coupling responsible for the non-linear effects 
comes from the long wavelength modes with k < 0.1/i/Mpc, 
so it is reasonable to assume that the non-linear effects can 
be described with a change of amplitude and slope, scaling 
linearly with ttiv, 



0.01 



0.15eV 



(23) 



The change of amplitude, crg_e£E = erg — 0.04(m;^/0.15eV), 
is already automatically included since we compute ag for 
a given cosmological model using its power spectrum. The 
spectral index seen by the non-linear correction is not the 
actual one, but the effective one given by the above equa- 
tion. This is justified by noting that the most important 
quantities that determine the shape of the non-linear cor- 
rection are the amplitude and slope of the power spectrum 
at a relevant pivot scale (fc ~ 0.1/i/Mpc in our case). Since 
the change in amplitude is already included in the change of 
erg, we just approximate the change in the slope of the linear 
power spectrum at the pivot point for m„ as the change in 
the spectral index. 

To test this procedure, we compare the resulting non- 
linear to linear power spectrum correction for massive neu- 
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trinos to the full simulations presented in iBird et alj (|2012h . 
For example, for the total neutrino mass of m,^ — 0.15 eV, 
we find that the reduced linear amplitude of the power spec- 
trum at the pivot point is 8 per cent, corresponding to a 4 
per cent change in ag. This leads to a further non- linear sup- 
pression of power, up to 3 per cent at fe ~ 1/i/Mpc. In addi- 
tion, the effective slope is reduced by 0.01 at the pivot point. 
This means that for massive neutrinos there is more power 
on large scales than in the zero mass case, relative to the 
pivot point. As a result, there is more mode-mode coupling 
which increases the small scale non-linear power, countering 
the effect from the reduced linear amplitude. The net effect 
is that the non-linear correction peaks at fc ~ 1/i/Mpc, but 
this quickly reverses sign above 2/i/Mpc. The overall effect is 
in a good agreement with the results of the full simulations 
of massive neutrinos in lBird et al. This suggests that 

we can parametrise the non-linear effect of massive neutri- 
nos simply with the change in the effective amplitude and 
slope of the linear power spectrum. 

To summarise the neutrino mass effects: at ~ 
O.S/i/Mpc, which is the peak of the contrib ution to T at 
R = 5fe~^Mpc and Ro = 3ft"^Mpc (Fig. 2 of iBaldauf et all 
I2OIOI ) and where we expect to have the most stringent con- 
straints from our data set, we expect about 10 per cent sup- 
pression of the power for "^m,^ = 0.15eV, relative to the 
zero mass case. 



2.3.3 Cosmology corrections 

When estimating the signal from the data, we assume a cos- 
mological model in order to convert angular distances A9, 
shears 7*, and redshift-space separations Az to transverse 
separation R, lensing surface density contrast AE, and line- 
of-sight separation H. Thus, for the model predictions for 
any other cosmology than the fiducial cosmology, we should 
in principle include a factor in both the transverse separa- 
tion and the amplitude of the measured signals to account 
for the fact that the wrong cosmology was used to do these 
conversions from observed to physical separations. However, 
for the highest redshift sample (for which this is most impor- 
tant), the size of the corrections is typically < 1 per cent for 
the range of allowed cosmological models. The correction is 
even smaller for the other samples, therefore it is well within 
the statistical error for this analysis, and we do not include 
it. 



2.3.4 Combining the model ingredients 

Finally, we need to combine these model ingredients to ob- 
tain ^gm(T') and Cgg(r). We do so by using the non-linear 
matter power spectrum for a given cosmology from Sec. 12.3.11 
along with Eos. [TSl and [T9l 

The results are then integrated to obtain the projected 
statistics that we use in reality. For the lensing signals, we 
integrate the correlation function along the line-of-sight to 
±140/i~^Mpc, consistent with the fact that the lensing win- 
dow is extremely broad. We do not include that window 
directly, but as shown in iBaldauf et all l|2010l ) fig. 9, its ef- 
fects are very small on the scales we use for science, and can 
be corrected for in a single factor that includes the cluster- 
ing line-of-sight integration length, redshift-space distortions 



(RSD), and the lensing window. Given that this correction 
factor is ~ 3 per cent at 60/i~^Mpc, much smaller than the 
observational errors, and ~ 1 per cent below 30/i~^Mpc, we 
neglect this correctior[f|. For clustering, we integrate along 
the line-of-sight to ±60/i~^Mpc, consistent with the observa- 
tional measurements. Finally, for both the lensing and clus- 
tering we use the projected surface densities to obtain T. 

Because we have some uncertainty in the calibra- 
tion of the lensing signal due to several systematic errors 
(Sec. I4.2~T1) . we include a nuisance calibration bias param- 
eter for the g-g lensing, which is assumed to have a mean 
zero and a Gaussian width of (4, 4, and 5) per cent for the 3 
samples. The calibration bias is assumed to be the same at 
all radii and for all samples - i.e., if the calibration bias is 4 
per cent for Main-L5 then it is 4 and 5 per cent for LRG and 
LRG-highz. This treatment is appropriate since the lensing 
calibration biases originate from the same sources for each 
sample, and improper estimation and removal of those biases 
would affect all samples nearly equally. 

2.4 Choice of Ro 

IBaldauf et ahl l|2010l ) considered relatively large values of _Ro 
such as 3 and 5/i~^Mpc. Use of a large value of Ro is advan- 
tageous from the perspective of systematic error, because it 
means that we are less sensitive to several effects that tend 
to be worse at small scales: cross-correlation coefficient de- 
viations from 1, deviations from our model for non-linear 
bias, and observational issues such as intrinsic alignments. 

However, use of larger _Ro will necessarily remove more 
of the measured signal, resulting in a noisier measurement. 
We therefore revisit the choice of Rq in order to achieve 
a fair compromise between statistical an d systematic error. 
Moreover, unlike in IBaldauf et al.l (|2010l ). we permit differ- 
ent Rf) for the galaxy-galaxy lensing and the galaxy cluster- 
ing, which is possible in the case that we explicitly model 
the signals (i.e., using different Ro means that results can- 
not be obtained by taking ratios of the two signals). The 
galaxy clustering signal is more sensitive than the galaxy- 
galaxy lensing to the fidelity of the non-linear bias model 
(because it has higher signal-to- noise) , so it might require a 
higher Ro to avoid systematic errors. 

We choose values of Ro for the two measurements based 
on modeling the simulated LRG sample using the same ma- 
chinery that we use to model the data (but without adding 
lensing shape noise, so that deviations from the real cos- 
mology in the simulations are due to a real analysis bias). 
We decided to use Ro = 2 and 4ft~'^Mpc for galaxy-galaxy 
lensing and galaxy clustering, respectively, since we find this 
choice still gives a small systematic error compared to the 
statistical error, as shown in Appendix |X] The correspond- 
ing plots of T are shown in Fig. [T] and we see that the 
simulated data and the best model agree reasonably well 
both for galaxy-galaxy lensing and galaxy clustering, with 

^ Technically, we have only done this test for the LRG sample. 
However, for the higher redshift and mass sample, the galaxy bias 
is higher and therefore RSD are even less important. For the lower 
redshift and mass sample, while the galaxy bias is lower and RSD 
are more important, we will see that the observational errors are 
also larger than for LRG. 
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62 ~ 0.25 as the best fit value. Note that the same value of 
62 also fits the cross-correlation coefficient for ^, as seen in 
the top panel of Fig. [1] 



2.5 Parameter values and constraints 

Our convention is to quote parameter values based on the 
median of the posterior distribution after marginalisation 
over all other parameters. The quoted error-bars come from 
using the PDF to determine the 68, 95, and 99.7 per cent 
confidence intervals. 



3 DATA 

Here we describe the data used in this paper, all of 
which come s from the Sloan Digital Sky Survey (SDSS). 
The SDSS l|York et all bOOOD imaged roughly tt steradi- 
ans of the sky, and followed up approxim ately one million 
of the detected objects spectroscopically ( Eisenstein et al] 
I2OOII : [Richards' et al. I I2OO2I : IStrauss et al.ll2002l ). The imag- 
ing was carried o ut by drift-scann i ng the sky in ph oto- 
metric conditions (iHogg et al. 200 ll: Ivezic et al. 2004), in 
five bands (ugriz) (|Fukugita et al.lll996l : ISmith et al."200j) 
using a specially-designed wide-field camera (Gurm c t al.l 
Il998l ). These imaging data were used to create the clus- 
ter and source catalogues that we use in this paper. All of 
the data were processed by completely automated pipelines 
that detect and measure photometric p roperties of objects . 



and astrometric ally calibrate t he da ta (|Lupton et al.l|200l| ; 
iPier et al.ll2003l : pTucker et all 1200^ . The SDSS-I/II imag- 
ing surveys were com pleted with a seventh data release 
ijAbazaiian et al]l2009l l. though this work will rely as well 
on an improved data reduction pipeline that was part of 
the eighth data release, from SDSS-III (Aihara et al. 2011); 
and an improved photome tric calibration ('ubercalibration', 
iPadmanabhan et al.|[2008l '). 

Below we describe the samples that are used as lenses 
and as sources. 



3.1 Main sample lenses 

The first lens sample that we use for this work is the 
flux- limited Main galaxy sample l|Strauss et al.l [20021 ) from 
SDSS DR7. The nominal flux limit is r < 17.77, when de- 



(based on a modifica- 



fined using Petrosian magnitude^ 

tion of | Petrosianlll976l described in Blanton et al.ll200ll and 
lYasuda et al.l 2001 ). In reality, the actual flux limit varies 
slightly across the survey area. We use the Main sample se- 
lection from t he NYU Value-Added Galaxy Catalog (VAGC, 
iBlanton et al.i .2005). which includes 7966 deg^ of spectro- 
scopic coverage (though we will employ further area cuts, 
described below). 

We select our sample using the 'LSS sample' DR7-2 
in the VAGC, which carefully tracks the spectroscopic flux 
limit and completeness across the sky. The particular LSS 



^ All magnitudes quoted in this paper are corrected f or Ga lactic 
extinction using the dust maps from 'Schlegel et al.l lll998l) an d 
the extinction- to-reddening ratios from lStouehton et al. 1 120021) . 



samples that we use are 'dr72full8' through 'dr72fulll0', 
where 'full' samples have the following properties: 

• They use all galaxies from r > 10 to the position- 
dependent flux limit. 

• They use areas with any level of completeness (even 
< 0.5, which occurs very rarely). 

• Galaxies that did not get a spectrum due to flbre- 
coUisions are assigned a redshift using the nearest-neighbour 
method. 

The '8', '9', and '10' subsamples have the following proper- 
ties, some of which will be subject to more cuts described 
below: 

• Redshift 0.001 < 2 < 0.4 

• The fc-corrections a re to z = 0.1 (kcorrect v4_1_4; 
iBlanton fc RoweidlioOTl ). 

• The distance modulus /j, — 5 log [Di/(10pc)] is deter- 
mined using Qrn ~ 0.3, flA ~ 0.7. This is formally inconsis- 
tent with the numbers used in the rest of this paper. How- 
ever, this is not a signiflcant issue here where we simply seek 
a reasonably volume-limited and consistent sample (partic- 
ularly given the weak dependence of the distance modulus 
on cosmology for these redshifts). 

• The luminosity evolution is assumed to have the form 

M{z) = M{z = 0.1) - 2[1 - (z - 0.1)](z - 0.1). (24) 

• The absolute magnitude is defined, in terms of the Pet- 
rosian magnitude r and galaxy redshift z, correcting for lu- 
minosity evolution, as 



M,. = r - [/i + ko.i{z) + M{z) - M{z = 0.1)]. 



(25) 



Given that the luminosity evolution is such that galax- 
ies were brighter in the past, the definition here removes 
that trend, connecting galaxies at one redshift to those 
that were suitably brighter at earlier times according to the 
empirically-determined evolution law in Eq. (|24p . 

• The absolute magnitude is then required to be in the 
range [-22,-21], [-21,-20], and [-20,-19] for the three 
samples, respectively. 

The effective area of the LSS sample is 7279 deg^ We 
then imposed some additional cuts on the area, removing 
regions without any source galaxi es in the backgr ound or 
within a Tycho bright star map l|H0g et al.ir2000l ). These 
cuts reduce the effective area to 7131 deg^. 

We wish to avoid overlap with the LRG lens samples 
described in the upcoming subsection, so that the cross- 
covariance between the signals with different lens samples 
can be assumed to be zero. Thus, we first require 0.02 < 
z < 0.155 (where the lower redshift cut removes galaxies for 
which it would be computationally prohibitive to measure 
correlations to 70/i~^Mpc, and the upper redshift cut re- 
moves overlaps with LRGs). W e then defined the three LS S 
samples using the notation from lMandelbaum et al.l (|2006al ): 
L3 with -19 > Mr ^ -20, L4 with -20 > Mr ^ -21, L5 
with —21 > Mr J? —22. We do not define any brighter sam- 
ples because their low abundance means that there are very 
few galaxies in those samples after the z < 0.155 cut. 

In practise, carrying out the analysis with all three sam- 
ples requires caution: if the redshift ranges overlap (as natu- 
rally occurs for a flux-limited sample) , then for scales above 
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Figure 2. Area coverage of the lens samples used in this paper. 



R ~ 8h~ Mpc, we find that the clustering and lensing sig- 
nals exhibit high cross-correlations between the samples - 
typically 80 per cent - because they trace similar large-scale 
structures. When limiting to volume-limited samples that do 
not overlap, the statistical power of L3 and L4 becomes rela- 
tively low on cosmological scales. In addition, our non-linear 
bias model in Sec. [2] was only tested on simulations with rel- 
atively high-mass halos (Mvir > 1 X 1 0" ft" ^ Ma). Given 
the typ ical halo masses for L3 and L4 in lMandelbaum et al.l 
l|2006al ). we conclude that the optimal way of including the 
Main sample in this analysis is to use L5 only. This sample 
(referred to as 'Main-L5' in the rest of the paper) includes 
69 150 galaxies, and is volume-limited for the redshift range 
that we use. 

For computation of cosmological observables, we require 
a set of random points that are distributed in the same way 
as the lens sample. We therefore use the sets of random 
points distributed for this sample in the VAGC LSS sample; 
this includes a proper redshift distribution that depends on 
the position-dependent flux limit at each point. For weight- 
ing, we use the inverse of the 'sector completeness' defin- 
ing the redshift success rate. The sector completeness for 
this sample has a median value of 0.972; 95 per cent of the 
galaxies have completeness above 0.924. 

The area coverage of the lens sample (7131 deg^, or 
/sky = 0.17) is shown in Fig. [51 This c overage is str i ctly a 
subsample of the source catalogue from iReves et all (|201ll ) 
that we describe in Sec. 13.31 

The redshift distribution dp/d2: and the comoving num- 
ber density n are shown for the Main L5 sample, and for the 
other samples described in subsequent subsections, in Fig. |3l 
For the Main-L5 sample, the comoving number density is 
~ 10-=*(/i/Mpc)^ 

The properties of this sample (and those introduced in 
subsequent sections) are summarised in Table [2l For the 
lensing, the effective redshift Zcff ,gm is determined not just by 
the lens redshift distribution, but also by geometric factors 
related to the relation between lens and source redshifts that 
come into the weighting scheme we use for estimating the 
signals (Sec. |4|. 



15 - 




z 



Figure 3. Top: Redshift distribution dp/dz for the three lens 
samples used in this work, as labeled on the plot. The solid 
and dashed lines show the unweighted and weighted histograms, 
respectively. Bottom: Comoving number density ft, in units of 
10-*(h/Mpc)3 (multiplied by 0.1 for Main-L5 for easier view- 
ing). The dotted vertical lines show the delineation between the 
different lens samples. 

3.2 LRG sample lenses 

We also define two lens samples consisting of Luminous Red 
Galaxies, or LRGs (Eisenstein et al.. ,2001i ). These galaxies 
have been used for numerous cosmology analyses with SDSS, 
most notably the detection of Baryon Acoustic Oscillations 
(BAO), which is enabled by the high galaxy bias of ~ 2 and 
the large volume probed by this sample, ~ 0.65{h~^Gpc)'^ . 

For selectio n of Luminou s Red Galaxies, we follow the 
methodologjQ of iKazin et al.l (|2010l ). which also starts from 
the NYU VAGC LSS sample described in the previous sub- 
section. In this case, only regions with completeness ^ 0.6 
are included; this definition is inconsistent with that used 
for Main-L5, but in practise, the discrepancy only affects 13 
deg^, or 0.2 per cent o f the area. Our sele ction is otherwise 
identical to that from iKazin et ahl (|2010l ). with the excep- 
tion of area cuts to restrict to the footprint of the source 
catalogue, eliminating 8 per cent of the LRGs. 

Rest-frame absolute magnitudes in the g band are cal- 
culated starting from the r band extinction-corrected ap- 
parent Petrosian magnitude. The distance modulus assumes 
f^m ~ 0.25 , SIa ~ 0.75. fc-c o rrecti ons and evolution correc- 
tions from lEisenstein et al.l HooJ) are used to convert r to 



iKazin et al] l|2010t ) have a well-defined procedure for 
calculating weights, completeness factors, dealing with fi- 
bre collisions, and distributing random points. In brief, they 
begin with a calculation of sector completeness that takes 
into account fibre collisions (i.e., if a galaxy did not get a 

http : //cosmo . nyu . edu/~ eak306/SDSS-LRG . html 
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Sample 




z range 




■^eff ,gm 


n 


name 










[ih/Mpcf] 


Main-L5 


69 150 


0.02 <: z < 0.155 


0.115 


0.109 


10-^ 


LRG 


62 081 


0.16 ^ 2 < 0.36 


0.278 


0.257 


10-4 


LRG-highz 


35 088 


0.36 ^ 2 < 0.47 


0.407 


0.398 


[8 - 6((2 - 0.36)/0.11)] X 10-5 



Table 2. Properties of the lens samples described in Sec. [3] For each sample, we show the redshift range, the mean redshifts, the effective 
redshift for the lensing as described in the text, and the number density. 



spectrum because it was not allocated a fibre, then this is 
not considered to be incompleteness). This completeness is 
used when distributing random points in the survey area; in 
any given region, they are diluted according to that region's 
sector completeness. To deal with the ~ 2 per cent of LRG 
targets that were not allocated fibres due to fibre collisions, 
a special weight is assigned; e.g., in a group of 3 LRGs of 
which only 2 were allocated a fibre, those two would each 
get a weight of 1.5. The random points - of which there are 
fifteen times as many as real points - are assigned a random 
redshift drawn from the p{z) of the real LRGs. 

We define two redshift samples, which we call 'LRG' 
(0.16 ^ z < 0.36) and 'LRG-highz' (0.36 ^ z < 0.47). In 
both cases, the absolute magnitude limits are —23.2 < Mg < 
—21.2; the former is approximately volume-limited, whereas 
the latter is flux- limited but relatively narrow (see Fig. |3]). 
In the flrst case, we define a radial weighting factor that is 
defined by the ratio of a smooth model of the p{z) to the 
real p{z). This weight helps eliminate the impact of minor 
fluctuations in the LRG density with redshift. Thus, for the 
LRG sample, the weights used for real points are 



fibre collision weight 



completeness 



Pmodcl(z) 



p{z) 



(26) 



and for random points, the same but without any fibre col- 
lision weig htfl For the LRG-highz sample, we use 



f'LRG-highz = 



fibre collision weight 
completeness 



(27) 



Once we include the redshift-dependent weighting, the 
LRG sample has a comoving number density that is nearly 
constant at 10~''(/i/Mpc)'^, a factor of ten smaller than for 
Main-L5. The LRG-highz sample has a comoving number 
density that drops with redshift because the sample is flux- 
limited. More details of these samples are shown in Tabled 



3.3 Sources 

The catalogue of source galaxies wi th measured shape s 
used in this paper is described in iReves et all (|201lD . 
hereafter Rll, which uses the re-Gaussianization method 
l|Hirata fc SeliakI 120031 ) of correcting for the effects of the 
point-spread function (PSF) on the observed galaxy shapes. 
The treatment of systematic errors is updated and improved 
compared to the previou s SPSS source catalogue u sing this 
PSF-correction method l|Mandelbaum et al.ll2005l ). in part 



° When normalising ratios of real versus random points, we use 
weights rather than absolute numbers of galaxies, and must watch 
out for the fact that if N^^al = ^rand. E^'roal Zl^'rand. be- 
cause of the fibre collision weighting on the real points. 



using tests of simulated SDSS images using real galaxies 
from COSMOS and real SDSS PSFs (|Mandelbaum et all 
'2013). To estimate source redshifts, we use photometric red- 
shifts (photo-z) based on the five-band photometry from 
the Z urich Extragalactic Ba yesian Redshift Analyzer (ZE- 
BRA, 'Feidmam|^^|200i), which were characterised by 
Nakaiima ct al.1 l|2012h . hereafter N12. 

The catalogue production procedure was described in 
detail in Rll, so we describe it only briefly here. Galax- 
ies were selected in a 9243 deg'^ region, with an average 
number density of 1.2 arcmin"^. The selection was based 
on cuts on the imaging quality, data reduction quality, 
galac tic extinction Ar < 2 defined using the dust maps 
from ISchlegel et al.l (Il998l) and the extinction-to-reddening 
ratios from IStoughton et al.1 l|2002D . apparent magnitude 
(extinction-corrected r < 21.8), photo-z and template used 
to estimate the photo-z, and galaxy size compared to the 
PSF. The apparent magnitude cut used model magnitude^, 
which are defined by fitting the galaxy profile in the r band 
to a Sersic profile with n = 1 (exponential) and n — 4 (de 
Vaucouleurs) , choosing the better of the two models, and 
then using that same rescaled profile to get magnitudes in 
all the bands. For comparing the galaxy size to that of the 
PSF, we use the resolution factor R2 which is defined using 
the trace of the moment matrix of the PSF Tp and of the 
observed (PSF-convolved) galaxy image Ti as 



(28) 



We require R2 > 1/3 in both r and i bands. 

The software pipeline used to create this catalogue ob- 
tains galaxy i mages in the r and i fi lters from the SDSS 
'atlas images ' (|Stoughton et al.ll2002l ). The basic principle 
of shear measurement using these images is to fit a Gaus- 
sian profile with elliptical isophotes to the image, and define 
the components of the ellipticity 



(e+,ex) = 



(b/a) 



1 + (b/a) 



■(cos 2(/!>, sin 2(j)), 



(29) 



where b/a is the axis ratio and (j) is the position angle of 
the major axis. The ellipticity is then an estimator for the 
shear. 



(7+,7x) = ^((e+,ex)), 



(30) 



where TZ « 0.87 is called the 'shear responsivity' and 
represents t he response of the ellipticity (Eq. to a 
small shear (|Kaiser et al ] I1995I : iBernstein fc Jarvisl I2OO2I '): 



^ http : / /www . sdss3 . org/drS/algorithms/ 
magnitudes . php#iiiag_model 
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Figure 4. Histograms of source galaxy properties, derived from 
a random subsample of 5 per cent of the catalogue after imposing 
all cuts (~ 2 X 10® galaxies). Top: Histogram of r-band extinction 
corrected model magnitude. Bottom: photo-2 histogram, and the 
inferred true dA''/d2 from N12. 



7^ : 



In the course of the re-Gaussianization PSF- 



correction method, corrections are apphed to account for 
non-Gaussianity of both the PSF and t he galaxy surface 
brightness profiles (jHirata fc Sell aki [20031 '). 

For this work, we do not use the entire source catalogue, 
only the portion overlapping the aforementioned lens sam- 
ples and around the edges. Fig. [4] shows histograms of the 
source galaxy r-band apparent magnitude and photo-2:. 



4 OBSERVATIONAL METHOD 

In this section, we describe how we use the galaxy catalogues 
from Sec. [3] to measure our two observables, the galaxy- 
galaxy lensing and the galaxy clustering. 



4.1 Galaxy-galaxy lensing 

Here we describe the computation of the lensing signal. 
For this computation, we rely on the lens catalogues in 



Sections 13.11 and 13.21 and the catalogues of random lenses 
with the corresponding area coverage and redshift distri- 
butions. First, pairs of lenses and sources that are physi- 
cally close on the sky and satisfy Zs > zi (using photo-z 
for sources) are identified. Here, "physically close" is deter- 
mined using the comoving transverse separation at the lens 
redshift; we require 0.02 < R < 32.9/i"^Mpc for Main-L5, 
and 0.02 < R < 73.1/i"^Mpc for the two LRG samples. 
These ranges are split into 37 or 41 logarithmic radial bins 
with Aii/i? = 0.2. 

Lens-source pairs are assigned weights according to the 
error on the source shape measurement via 



Wis 



siK + a|^) 



(31) 



where cr|jv, the intrinsic shape noise, was found in Rll to be 
0.365, independent of magnitude. The factor of E^^ means 
that we weight by inverse variance of the expected lensing 
signal, thus downweighting pairs that are close in redshift 
because the lensing geometry is suboptimal. 

Once we have computed these weights, the lensing sig- 
nal in each annular bin can be computed via a summa- 
tion over lens-source pairs "Is" and random lens-source pairs 
"rs": 



AE(i?) = 



(32) 



where the factor of 27?. arises due to our definition of elliptic- 
ity. Note that t his is equivalent to the pro cedure in previous 
works such as iMandelbaum et al.l (|2005h of using 'Yl.is 
in the denominator and then multiplying the result by the 
boost factor. 



B{R) = 



J2 ^rs 



(33) 



The division by Wrs is necessary to account for the fact 
that some of our 'sources' are physically associated with the 
lens, and therefore not lensed by ilF°l. 

Due to the observational strategy in SDSS, there is a 
tendency for the PSF to align coherently along the scan 
direction. This tendency gives rise to a so-called 'system- 
atic shear' if the PSF correction is not perfectly efficient at 
removing the PSF ellipticity from the galaxy shapes, which 
turns out to be the case at a low level for our PSF-correction 
method. If a lens has a uniform distribution of sources 
around it, the contribution of the coherently- aligned sys- 
tematic shear to the average tangential shear is zero. Thus, 
the systematic shear can contribute to the lensing signal 
primarily due to the inclusion of lenses near survey bound- 
aries, since they lack a symmetric distribution of sources. 
To remove this systematic shear, we can simply subtract 
the lensing signal AE^and measured around random points, 
which will capture the geometry-dependent effect of th e sys- 
tematic shear. As noted bv iMandelbaum et al.l l|2005h . this 
correction may be imperfect in the case that the lens den- 
sity fiuctuates due to some effect that also modulates the 



This correction is formally correct in the limit that other ef- 
fects that modulate the source number density, such as magnifi- 
cation or difficulty in detecting sources due to software or light 
from the lens, are negligible. This condition is satisfied on all 
scales used in this work. 
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systematic shear, if this modulation is not included in the 
random point distribution. We will return to this issue in 
Sec.lH 

To compute T using Eq. (|12[l on our noisy, binned data, 
we use the following procedure. We first determine AE(i?o); 
a discussion of how this procedure can affect results is in 
iMandelbaum et al.l (j2010l ). For this project, as we will show 
in Sec. IBll we are helped by the fact that around Ro there 
is a range of scales on which AE(J?) is well-approximated 
by a power law. This si tuation is different from that of 
IMandelbaum et al.l (|2010l ). which used galaxy clusters with 
Ro well within the cluster virial radius, so that the scaling 
of AS(_R) was inconsistent with a power law. Thus, in this 
paper (unlike previous work) we can fit AE to a power law 
over a fixed range of scales, which will minimise the noise 
in estimation of AE(i?o) that gets propagated into T. We 
present details and tests of this procedure in Sec. IBll After 
estimating AE(_Ro), we compute T in each radial bin using 
Eq. directly. 

For science, we only use scales above 2/i~^Mpc 
(4/i^^Mpc) for the lensing (clustering), which results in the 
inclusion of 18 radial bins for LRGs and 14 for Main sam- 
ple lenses (14 and 10). These bin counts take into account 
the fact that we exclude the nominal first radial bin above 
4/i~^Mpc for the clustering, because while the bin center is 
above _Ro, the lower edge is below Ro, and we do not at- 
tempt to model T in this rapidly-varying regime. The max- 
imum scales of ~ 30 and ~ 70/i~^Mpc for Main-L5 and the 
LRG samples, respectively, were chosen based on consider- 
ations related to systematic error, which will be described 
in Appendix IB5I In brief, our finding is that there seem to 
be fluctuations of lens number densities that correlate with 
systematic errors in the shear for larger scales and lead to 
a situation where our systematic uncertainty exceeds the 
statistical error on the signal. 

To determine errors on the lensing signal AE or derived 
quantities like T, we divide the survey area into 100 equal- 
are j"! iackknife subregions, each of size ~ 71 deg^ or typical 
length scale ~ 8.4 degrees. The same division of the area into 
regions will be used when computing the galaxy clustering 
signal, so that we can also estimate the covariance between 
the two. This number of regions was motivated by a desire 
to balance two competing effects. First, we require that the 
number of regions be significantly larger than the number of 
radi al bins (18), to red uce the noise in the covariance ma- 
trix l|Hirata et al.|[20o3 ). Second, we require that the region 
size be larger than the maximum scale used for science. For 
the Main sample, 30fe~^Mpc (comoving) at z — 0.11 cor- 
responds to 5.3 degrees; for LRGs, 70/i"^Mpc at z = 0.27 
corresponds to 5.2 degrees. Thus, our typical region size is 
60 per cent larger than the maximum angular scale used for 
science. 

When computing the covariance for derived quantities 
such as T, we estimate T for each jackknife sample to get the 
covariance matrix, rather than using the covariance matrix 
for AE and propagating errors. 

It is well known that jackknife covariance matrices can- 



For an arbitrary survey geometry, it is difficult to achieve 
equal-area and contiguous regions. We have opted for equal-area 
regions, roughly 10 per cent of which are not contiguous. 



not be used to get cosmological constraints without some 
correctio n due to the fini te level of noise (e.g.. iHirata et all 
i2004: Hartlap et al.ll2007l l: this is a consequence of the fact 
that the inverse of a noisy, unbiased estimator of the co- 
variance matrix is not an unbiased estimator of the inverse 
covariance matrix. We handle this issue by modeling the co- 
variance matrix to eliminate noise. (While this might seem 
to eliminate the need to make many jackknife regions to 
reduce noise, as we have already done, we still need the co- 
variance matrix to be reasonably well-determined in order to 
easily model it empirically.) Details of this approach will be 
described in Sections l5.2l and l6.2l However, we note that our 
results are insensitive to whether we use the noisy jackk nife 
covariances with a correction factor (jHartlap et al ][2003) af- 
ter inverting to obtain the inverse covariance, or whether we 
use the covariance matrices that we have modeled to reduce 
the noise. This finding suggests that our results are not sig- 
nificantly impacted by systematics related to our handling 
of covariance matrices. 



4.2 Lensing systematic errors 

A thorough treatment of systematic errors with this source 
catalogue is in Rll. Here we include only a brief summary 
of the issues, along with the impact for this work. 

4-. 2.1 Calibration biases 

In iReves et al.1 (|201ll ) we considered several different types 
of systematic errors for which we applied corrections and 
estimated a total error budget. In this work we consider the 
same set of systematic errors, with the only change being 
that the lens samples are at different redshifts, thus chang- 
ing the values of many of the systematic errors and their 
uncertainties. 

To summarise briefiy, our approach to estimating the 
systematic error budget is to consider a full list of systematic 
errors that affect the lensing signal calibration. We correct 
for our best estimate of any biases, and assign systematic 
errors using the following prescription: for those types of 
biases that are inherently connected, we assume that sys- 
tematic uncertainties add linearly (e.g., two sources of 1 per 
cent-level uncertainty become a combined 2 per cent un- 
certainty); for those that are independent, we add them in 
quadrature (i.e. in the previous example, the combined un- 
certainty would be \/2 per cent). 

There are three calibration biases related to shear esti- 
mation that we consider to be inherently connected: errors 
in the correction for PSF dilution due to the PSF correction 
method not being perfect; noise rectification bias; and selec- 
tion biases (due to our resolution cut favouring galaxies that 
are aligned with the shear). In Rll we described tests us- 
ing realistic galaxy simulations ijMandelbaum et al.ll201 to 
constrain these three uncertainties together, which yielded 
a combined 3.5 per cent uncertainty largely independent of 
the lens redshift. 

The other calibration biases that we consider to be inde- 
pendent are the impact of photo-z error (as discussed thor- 
oughly in N12); stellar contamination, which we constrain 
using space-based data; PSF model uncertainty; and shear 
responsivity errors due to incorrect estimation of the RMS 
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galaxy ellipticity. Of these, the first is the dominant one (1, 
2, and 3 per cent uncertainty for Main-L5, LRG, and LRG- 
highz respectively - because as shown in N12, the systematic 
uncertainty is larger for higher redshift samples, where cos- 
mic variance in the calibration samples is more important). 
Both stellar contamination and PSF model uncertainties are 
< 0.5 per cent. Shear responsivity uncertainty is 1 per cent 
for all samples. Thus, the three shear biases listed previously 
are the dominant uncertainty for all samples; when we add 
up the independent effects in quadrature, we obtain a 4, 5, 
and 5 per cent systematic uncertainty for Main-L5, LRG, 
and LRG-highz, respectively. 

For the purpose of simplifying the modeling, we assume 
that this final calibration uncertainty has a Gaussian error 
distribution, which may not be quite correct in detail. More- 
over, since the errors were assessed in the same way for each 
lens sample, we assume that they are 100 per cent correlated 
- i.e. if the calibration is really 4 per cent too high for Main- 
L5, then it is 5 per cent too high for LRG and LRG-highz. 
We include this calibration uncertainty in the modeling of 
the lensing signal. 

To test our understandi ng of the calibration biase s, we 
present several ratio tests l|Mandelbaum et al] 120051 '). i.e. 
comparisons of the signal computed using the same lens sam- 
ples, but with different subsamples of the source catalogue. 
After we correct for our understanding of the calibration 
biases, we should find that the ratios of these signals are 
consistent with 1 within the errorcj. 



4-2.2 Scale- dependent systematics 

iMandelbaum et al.l (|2010l ) includes a list of scale-dependent 
systematic errors that complicate the inference of cluster 
masses from the cluster lensing signal. Fortunately, many 
such errors are sufficiently small for galaxy-scale lenses 
and/or on the > 2/i~^Mpc scales that we use for science that 
we can ignore them. The scale-dependent systematic errors 
that we do consider are intr insic alignments of galaxy shapes 
(e.g.. lHirata fc Seliaklliooih . given that we know some of our 
'sources' are really physically associated with the lens and 
therefore may tend to point towards the lens. In principle, 
this effect can be quite important if we have no way of re- 
moving galaxies that are physically associated with lenses 
from our source sample; fortunately, our photo-0 are suf- 
ficiently good that we are fairly successful at doing so. In 
Sec. IB21 we will estimate the importance of this effect based 
on the fraction of ph ysicall y-associated galaxies as a function 
of scale (see also iBlazek et al.. 2012 ). 

The other main scale-dependent systematic error is the 
'systematic shear' described in Sec. 14. II "While we can use the 
procedure described there to correct for it, we also must test 
the validity of that correction procedure, which we will do 
once we present the results in Sec. lB5l Moreover, the system- 
atic shear is the main factor that determines the maximum 
scale that we use; our maximum scales of ~ 30 (Main-L5) 



and ~ 70/!.~^Mpc (LRG, LRG-highz) are motivated by a de- 
sire to avoid a situation where the correction for systematic 
shear is comparable in size to the real lensing shear. 

4.3 Galaxy clustering 

We compute the galaxy clustering signals using the same 
logarithmic binning size and maximum R as for the lensing, 
but with a minimum R — 0.3^~^Mpc (which provides some 
measurements below i?o for estimating Wgg(_Ro = 4)). 

The estimation of clustering signals for the lens samples 
relies on SDSSpijQ software to rapidly identify galaxy pairs 
within the required separation on the sky. To compute the 
galaxy auto-correlation w;gg(i?), we begin by computing the 
3D galaxy correlation function ^gg on a grid of values in 
{R, n) where 11 is the comoving line-of-sight separation with 
respect to the mean position of the galaxies in the pair. Our 
estimator for the correlation fun ction is a generalisation of 
that from lLandv fc Szalavl (|l993l ). 

t TT^ DD~2DR + RR 

Cgg (-R, n) = — , (34) 

using sums of products of weights rather than numbers 
of pairs of data-data, data-random, and random-random 
pairf^. Here, the weights for a given pair come from the 
product of the weight for each galaxy in the pair, where 
the weight per galaxy is initially defined as in Sec. [3] (e.g., 
Eq. I26[) . For the LRG and LRG-highz samples, there is an 
additional factor in the weight, to account for the fact that 
the g-g lensing and galaxy clustering measurements would 
have different effective weights since the g-g lensing auto- 
matically includes a lensing weight factor that depends on 
the redshift distribution of the source galaxies. This lensing 
weight is a decreasing function of redshift, and we include 
it in the clustering analysis so that the two measurements 
will not have different effective amplitudes (or even shapes, 
since the full scale-dependent matter clustering and non- 
linear bias can evolve with redshift). In practise, it turns 
out that this weight is not important for Main-L5 because 
those galaxies are well below the bulk of the source redshift 
distribution and because the Main-L5 redshift distribution 
is rather narrow. For LRGs, it is more important, chang- 
ing the effective redshift of the clustering measurement by 
Az = 0.03. 

The projected correlation function 'Wgg{R) is formally 
defined as 



(35) 



integrated along the entire line-of-sight. In practise, we com- 
pute Wgg{R) via a limited summation. 



,(7?) = An^Cgg(ii,n. 



(36) 



in 40 bins in 11 that are linearly spaced with AH = 3, span- 
ning a range -Ilmax ^ n < Ilmax, for Ilmax = 60/i~^Mpc 
(we consider the impact of this choice of Ilmax in Sec. [C}. 



For the combinations of lens and source redshifts used liere, tlie 
predicted differences in those ratios due to reasonable variations 
on our adopted cosmological model are at the 0.1 per cent level, 
well within the errors. 



http : //dls .physics .ucdavis . edu/~ scranton/SDSSPix/ 

To reduce the noise, we have many more random points than 
real points. Here, all numbers such as DR and RR (or their gen- 
eralisation in terms of pairs of products of weights) are properly 
normalised to account for this fact. 
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To estimate Tgg(_R), we use Eq. (|13|) . replacing E with 
Wgg, which requires an estimate of uigg(i?o)- As for the lens- 
ing signal, we identify a range of scales for which the signal 
is approximately a power law, and fit Wgg to a power-law on 
those scales to estimate Wgg{Ro). Tests of this determina- 
tion of Wgg(7?o) are presented in Sec. ICll We then employ 
Eq. (|13p to get Tgg; the first term, an integral over all scales 
above _Ro, is done via summation. For each bin, we effec- 
tively assume a constant w within the bin, and we carefully 
account for partial bins that fall in the R range of interest. 

As for the lensing signal, we use the division of the 
lens samples into equal-area jackknife regions to compute 
covariance matrices, and we present the results of the jack- 
knife covariances and modeling them to reduce the noise in 
Sec. [621 



4.4 Galaxy clustering systematic errors 

There are several possible systematic errors in the calcula- 
tion of the galaxy clustering statistics. 

One issue is the handling of the integral in Eq. (|36|) . with 
the finite line-of-sight cutoff. Naively, we can account for 
this in our modeling of the clustering signal by integrating 
the theoretical 3D correlation function to the same line-of- 
sight cutoff. However, as illustrated bv IPadmanabhan et al.l 
l)2007h . this approach is uncertain at the level of tens of per 
cent on large scales by redshift-s pace distortions (RSD), and 
a linear treatment ('Kaiser '1987") is not likely to be adequate 
on these scal es (e.g ., Rcid & White 2011). Fortunately, as 
illustrated bv lBaldauf et al.l (|2010l ). the uncertainty induced 
by the redshift-space distortions is far less important for the 
observable that we use for science, Tgg . As shown in figure 8 
(right panel) in that paper, the impact of RSD and the finite 
line-of-sight cutoff is reduced from ~ 30 per cent bias on uigg 
to ~ 5 per cent bias on Tgg at the maximum scale that we 
use for science, and is very small below 30/i~^Mpc. We could 
apply a correction for this small, residual systematic error, 
but as described in Sec. 12.3.41 a combined correction for this 
effect and others (e.g., the lensing window, which goes in 
the opposite direction) is so small as to be negligible for our 
analysis. 

Because Tgg is a partially compensated statistic, it is 
also less sensitive to large-scale density fluctuations that can 
shift Wgg up and down. We will demonstrate the effect of this 
on the cosmic variance component of the errors in Sec. 16.21 
The same argument is true for the integral constraint, which 
will lead to a constant offset in Wgg which goes away when 
computing Tgg. 

There are also a question of ho w sample de fi nition 
choices affect the measured statistics. iKazin et all (l2010l ) 
consider several such effects for uigg, including the method 
of distributing the random points in redshift, and the way 
of handling fibre collisions in the weighting. Their conclu- 
sion that these issues are only important at the few per cent 
level at most is also applicable to our results, and thus this 
systematic error is subdominant to the statistical and sys- 
tematic uncertainties in the lensing signal. 



^'lrg [degrees] 
0.01 0.1 1 




R [h-'Mpc] 



Figure 5. Top: Observed lensing signal ASgni(i?) for all three 
lens samples. The vertical line indicates the minimum scale used 
for cosmological constraints. The axes on the top and right indi- 
cate the angular scale 6 and the tangential shear 1007t for the 
LRG sample. Bottom: RTgrniR', Ro = 2) for all three samples as 
labeled on the plot. 

5 RESULTS OF LENSING MEASUREMENTS 

In this section we present the galaxy-galaxy lensing mea- 
surements (Sec. 15. ip . the error estimates (Sec. 15.2(1 . Tests 
for systematic errors are in Appendix [B] 

5.1 Observations 

The lensing signals for all three samples are shown in Fig. [5] 
This figure shows the observable quantity, AE, and also 
the quantity used for cosmological constraints, Tgm(^2o = 
2/i~^Mpc), plotted as RT for easier viewing on a linear scale. 
Clearly the S/N of the observable is quite high - typically 
~ 25 averaged over all scales, using 

|=pC-.)''^ (37) 

where x is the vector of AS values in each radial bin and 
C is their covariance matrix, to account for correlations be- 
tween bins. The shear is well-constrained down to a level of 
~ 5 X 10^^ (at ~ 5 degree angular separations), and the re- 
sults are statistically consistent with previously-pu blished 
ones f or LRGs with R < 4/i~^Mpc ( Mandclba um et all 
l2006bl ). Comparison with previous Main-L5 observations is 
complicated by our imposition of a redshift cut z < 0.155, 
and this is the first such galaxy-galaxy lensing observation 
for LRG-highz. 

Tgm is a lower S/N quantity, with total average S/N of 
11, 14, and 8 for Main-L5, LRG, and LRG-highz, averaged 
over the range of scales (2 < i? < 70/i~^Mpc) shown on the 
plot. 

A detailed discussion of tests for systematic errors in 
the lensing signal is in Appendix IbI 



© 0000 RAS, MNRAS 000, 000-000 



16 Mandelbaum et al. 



5.2 Covariance matrices 

Here, we present the error estimates for the Tgm results 
shown above. As stated in Sec. 14. II the noisiness of the jack- 
knife covariance matrices requires some correction in order 
to get cosmological parameter estimates. Rather than mod- 
ifying the procedure for using the m to get confidence i nter- 
vals, as in iHirata et al.l \2QQ4 ) or iHartlap et al.l ()2007l '). we 
instead model the matrices to make noiseless versions. 

The process begins by modeling the diagonal terms of 
the covariance matrix as a function of R. We refer to the co- 
variance matrix for Tgm as Cgm* with elements correspond- 
ing to radial bins i and j of Cgm {Ri, Rj)- Our smooth model 
is 

C^^J{R^,R^) = AR-\l + [R/Rtfl (38) 

with an amplitude A and a turnover radius Rt. This two- 
parameter model is motivated as follows: on all scales, we 
expect sampling variance to be minimal because of the large 
area and the compensated nature of T, so shape noise should 
be the dominant source of error. The shape noise variance 
scales like R~^ for logarithmically spaced annular bins, and 
we fit for the amplitu de A of this term. However, as shown 
in Ijeong et al.l (|2009l ). above some radius the shape noise 
fails to decrease as rapidly with R, in the regime where R 
is significantly larger than the typical separation between 
lenses. In that case, the lens-source pairs in the annular bin 
include many of the same sources around nearby lenses, so 
the shape noise does not decrease by adding more lenses. 
The term in brackets in Eq. (|38|) represents this flattening 
of the errors with scale. (There will also be a corresponding 
increase in bin-to-bin correlations on those scales, as will 
shortly be apparent.) 

Our approach is to model the diagonal elements of the 
covariance matrix by directly fitting for [A, Rt) for each lens 
sample using minimisation, doing an unweighted fit for 
log Cgm^(J?i, i?i) as a function of log J?. The scale on which 
the term in brackets in Eq. (|38[) becomes important is Rt — 
(8, 31, 41)ft~'^Mpc for the three samples. This trend of Rt 
is unsurprising given the trends in lens number density for 
the three samples. In the top panel of Fig. |S1 we show a 
comparison between the jackknife covariance diagonal terms 
and those from the model, for the LRG sample. As shown, 
the level of fluctuations of the jackknife variances compared 
to those in the model is 12 per cent. 

Next, we model the off-diagonal terms, which are also 
somewhat noisy. Off-diagonal terms can arise due to (a) 
cosmic variance (not very significant for this sample), (b) 
correlated shape noise due to the large R compared to the 
separation between lenses, and (c) the fact that T(_R) at 
some radius depends on the AE(_Ro), which tends to anti- 
correlate bins at 7? ~ _Ro with each other. Since there are sev- 
eral sources of correlations, they are not as simple to model 
analytically. Thus, we take a non-parametric approach, by 
generating the correlation matrix, i.e. Ccorr, defined by 

/-~c _ CgrajRi, Rj) /„(,N 

^corr,z,j — J V / 

Y Cgni {Ri, Ri ) Cgni' (-Rj , Rj ) 

We then apply a boxcar smoothing algorithm with a length 
of 3 bins in radius to this matrix, to reduce the noise. The 
middle and bottom panels of Fig. [6] show the unsmoothed 



and smoothed correlation matrix for the LRG sample. As 
shown, the smoothing has not resulted in any significant 
modification of the apparent real trends, but has eliminated 
the majority of the noise. 



6 RESULTS OF CLUSTERING 
MEASUREMENTS 

In this section we present the galaxy clustering measure- 
ments (Sec. [GTTJ, the error estimates (Sec. |6]2]), and cross- 
covariance with the lensing results (Sec. 16. 3) . Tests for sys- 
tematic errors in the clustering measurements are in Ap- 
pendix [C] 



6.1 Observations 

The clustering signals for all three samples are shown in 
Fig. [7] This figure shows the observable quantity, uigg, 
and also the quantity used for cosmological constraints, 
Tgg(i?o = 4/i"^Mpc) (plotted as RT for easier viewmg on 
a linear scale). Clearly the S/N of the observable is quite 
high, significantly more so than for the lensing observable. 
T gives a total average 5'/A'' of 19, 33, and 32 (Main, LRG, 
LRG-highz) when averaged over scales R > 4h~^Mpc using 
Eq. §7}. 

As discussed in Sec. 14.31 these results include a redshift- 
dependent weighting factor so that the effective redshift will 
be the same as for the galaxy-galaxy lensing measurement. 
Thus, they cannot be directly compared with previous mea- 
surements of LRG and LRG-highz sample clustering without 
checking the effect of this weighting. We find that for the 
LRG (LRG-highz) sample, inclusion of the lens- weighting 
factor has lowered the amplitude of Wgg by 4 (2) per cenll^^l 
and increased the errors by ~ 10 per cent. Thus, if we had 
not included it, then with the same lensing signal but a 
higher clustering signal we would have inferred a lower erg 
by 2 (1) per cent when analysing these samples, which is of 
similar order as other calibration factors we have considered 
and therefore validates our choice to include this redshift- 
weighting properly. 

The bottom panel of Fig. [7] shows the ratio of the mea - 
sured signal for the LRG sample from IZehavi et al.l (|2005l ) 
to our measurement without the lens-weighting, with the dif- 
ferences being due to our use of nearly twice as much area 
(7131 versus 3836 degrees^) and the different line-of-sight 
integration lengths (±60 versus ±80 h-^Mpc). Even with 
these differences, the results agree at the 1 per cent level. 
Additional tests for systematic errors in the clustering mea- 
surements are in Appendix [Cl 



The expected sign and magnitude of the effect is not com- 
pletely clear; for a passively evolving population the sign should 
in fact be the opposite. However, we have split the LRG samples 
into redshift slices and confirmed that within the redshift range of 
the LRG sample, the amplitude of Wgg evolves by as much as 10 
per cent, with lower amplitude at lower redshift, consistent with 
our findings with lens-weighting included. 
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Figure 6. Top: A comparison of the jackknife covariance matrix 
diagonal terms for Tgm for the LRG lens sample with the model 
covariance matrix terms as a function of transverse separation 
R. Middle: Jackknife correlation matrix for Tgm, again for LRG 
lenses. Bottom: Smoothed correlation matrix for Tgm. 



Figure 7. Top: Observed clustering signal uigg(R) for all three 
lens samples, including lens-weighting factors for LRG and LRG- 
highz. The vertical line indicates the minimum scale used for cos- 
mological constraints. Line colours and types indicate the sample, 
using the same scheme as Fig. |5] Middle: RTgg{R; Ro = 4) for 
all three samples as labeled o n the plot. B o ttom: For the LRG 
sample, the ratio of Wgg from lZehavi et aD l[200^ to that from 
our work without lens- weighting, after accounting for the different 
radial binning. 



6.2 Covariance matrix 

Here, we present the error estimates for the quantities shown 
above, in particular, the Tgg resuhs. We follow the approach 
from Sec. l5.2l of modeling the covariance matrices to reduce 
the noise, again showing results for the LRG sample as an 
example. 

The process begins by modeling the diagonal terms of 
the covariance matrix as a function of R. We refer to the co- 
variance matrix for Tgg as Cgg' with elements correspond- 
ing to radial bins i and j of Cgg''(7?i, Rj). Our smooth four- 
parameter model is 

Cg\R„R,) = AiR"' +A2R''\ (40) 

a sum of a shallower and a steeper power-law, which dom- 
inate on larger and smaller scales respectively. The choice 
of this functional form comes from the fact that for Tgg 
there is an additional noise component on small scales due 
to propagated uncertainty in ■u;gg(i?o)- 

We model this term by directly fitting for 
(j4i, ai, j42, Q2) for each lens sample using minimi- 
sation, doing an unweighted fit for log Cgg''(i?i, as a 
function of logT?^. In the top panel of Fig. |8] we show a 
comparison between the jackknife covariance diagonal terms 
and those from the model, for the LRG sample. As shown, 
the level of fluctuations of the observed variances compared 
to those in the model is 17 per cent (but there is a fairly 
extreme outlier; excluding that one noticeably decreases the 
estimated scatter). Smoothing the diagonal terms to reduce 
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the influence of this outHer is important; too low a variance 
would result in that bin unfairly dominating the fits. 

Next, we model the off-diagonal terms, which are also 
somewhat noisy. Off-diagonal terms can arise due to cosmic 
variance (not very significant for this compensated statis- 
tic Tgg), but also a small contribution from the fact that 
Tgg(_R) at some radius depends on Wgg(_Ro), which tends to 
anti-correlate bins at J? ~ _Ro with each other. Since there 
are several sources of these correlations, they are not as sim- 
ple to model analytically. Thus, we take a non-parametric 
approach, by generating the correlation matrix, i.e. Ccorr, 
defined by 



cg\R,,R,) 



(41) 



We then apply a boxcar smoothing algorithm with a length 
of 3 bins in radius to this matrix, to reduce the noise. The 
middle and bottom panels of Fig. [S] show the unsmoothed 
and smoothed correlation matrix for the LRG sample. As 
shown, the smoothing has not resulted in any significant 
modification of the apparent real trends, but has eliminated 
the majority of the noise. 

We also compare the covariance matrices for w^g and 
Tgg, to check that they behave in the way that we expect 
with respect to reduced cosmic variance in the latter due 
to its compensated nature (Sec. [2]). Fig. [9] illustrates this 
difference, using the LRG sample as an example. Here we 
have rebinned the data and covariances by a factor of two, 
because we do not have a smoothed version of the covari- 
ance matrix for Wgg. As shown, the relative error (top panel) 
on Tgg is larger near Ro, because Tgg is defined to be near 
zero there, but on large scales, uigg has a larger fractional 
error due to cosmic variance. We also study the correla- 
tion properties of these statistics. As shown (bottom panel), 
Wgg exhibits larger correlations between nearby bins. This 
is a consequence of cosmic variance: large-scale modes can 
coherently shift uigg up or down, resulting in bin-to-bin cor- 
relations. In contrast, Tgg shows less significant correlation 
patterns. 

We also compare these error properties versus the ex- 
pected ones from the simulations. For the simulations, the 
volume has been divided into 1.5 x 1.5 x 0.3 = 0.67(Gpc//i)^ 
sub-volumes, almost exactly the size of our observed vol- 
ume but with slightly different geometry. There are 40 such 
sub-volumes, which we use to estimate covariance matrices 
by comparing the signals computed in each one. This test 
serves as a check on the jackknife method that we carry out 
on the real data. We find that the scaling of cr{'Wgg)/'Wgg in 
the simulations and data is extremely similar for all scales 
that we use. For Tgg the same is generally true, though the 
errors near Ro seem to be ~ 20 per cent larger in the simu- 
lations than in the data, but the opposite is true at higher 
R, with the two in agreement around 8/i~^Mpc. This may 
be attributed to the method of determining Wgg{Ro), which 
differs slightly between the data analysis and simulations. 
Nonetheless, the comparison between errors in the simula- 
tions and in the real data is very similar, and should alleviate 
any concerns about the jackknife errorbars used in practise. 
We have also confirmed that the correlation properties for 

and Tgg in the simulations are consistent with those we 
observe in the real data (bottom of Fig. [9]). 
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Figure 8. Top: A comparison of the jackknife covariance matrix 
diagonal terms for Tgg for the LRG sample, with the model co- 
variance matrix terms as a function of transverse separation R. 
Middle: Jackknife correlation matrix for Tgg for the LRG lens 
sample. Bottom: Smoothed correlation matrix for Tgg. 
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Figure 9. Error comparison for Wgg (black solid lines) and Tgg 
(red dashed lines) for the LRG lens sample. Top: Standard de- 
viation of these statistics (i.e., the square root of the diagonal 
covariance matrix elements) as a function of R, divided by the 
signal itself. Bottom: Correlation of the data at each bin in ra- 
dius with the bin above R = i?o/i~^Mpc. 

6.3 Cross-covariance with lensing 

Using the jackknife resampling, we are able to compute the 
cross-covariance between the lensing and clustering observa- 
tions. Our findings are that on all scales, these correlations 
are consistent with zero, except possibly on the largest scales 
where they reach 10-15 per cent. This result is unsurprising 
given that the dominant source of error in our observations 
is the lensing shape noise. So, for the purpose of cosmolog- 
ical parameter constraints, we set the cross-terms between 
lensing and clustering in the covariance matrix to precisely 
zero. 

For future work with surveys that are not as limited by 
shape noise, these covariances will be important to model 
accurately. 



7 COSMOLOGICAL PARAMETER 
CONSTRAINTS 

In this section, we present cosmological parameter con- 
straints derived from the data that were shown in the pre- 
vious section. 

7.1 Constraints with these data alone 

We use the standard cosmomc (ILewis fc Bridlell2002l ) pack- 
age to make statistical inferences from our data. The package 
has been extended with a module that models our signals 
as described above and provides the likelihood that can be 
turned into posterior probability either alone or in conjunc- 
tion with another dataset. We use the standard cosmomc 



parametrisation of cosmology with flat priors on the follow- 
ing parameters ujt, — flt,h^ (baryon density), ujcdm (cold dark 
matter density), 9 (angular diameter distance of the surface 
of last scattering, a proxy for Bubble's constant), r (optical 
depth to the surface of last scattering), w (dark energy equa- 
tion of state). Us (spectral index of primordial fluctuations), 
A (amplitude of the primordial fluctuations; however, we 
quote our results in terms of the more commonly-used pa- 
rameter (Tg). In addition we vary parameters specific to our 
dataset: two bias parameters for each dataset and a common 
calibration parameter c. 

We start by performing the analysis using lensing data 
alone, both in individual redshift bins and combined. For 
these tests, we fix the cosmological parameters to ilm = 
0.25, ris = 0.96, and h — 0.7, only varying erg, the bias 
parameters, and the lensing calibration parameter. We do 
this fit for each sample separately, and for all three together 
(Table [1 fits 1-4). 

With these fits, we can do a basic sanity test for con- 
sistency between samples. In Fig. [10] we show the data with 
the best-fitting signals. The model appears to fit the data 
quite well, without any signs of systematic tension, consis- 
tent with the reasonable and p- values in Table [S] 

Fig. [11] shows the posterior distributions for the bias 
parameters for each sample; we show the results from fits 
1-3 (fitting for the samples separately) and 4 (fitting all to- 
gether). Fig. 1121 shows the ID posterior distributions for erg 
for each sample separately and the combined result (fits 1- 
4), in all cases marginalised over the bias parameters. The 
corresponding limits on erg when using all samples together 
(as in Table [1 fit 4) are 0.80 ± 0.05 (68 per cent confi- 
dence level, stat. + sys.). Figures [TCTlll II and ll2l demonstrate 
clearly that there is no statistically significant tension be- 
tween cosmological parameter constraints from the different 
samples. 

The results in Fig. [TT] suggest that for LRG and LRG- 
highz, we detect non-linear bias at the > 3cr and ~ 1.5a 
levels, respectively (the results in Table [3] are deceptive; 
they suggest less significant detections, because the error- 
bars are quite non- Gaussian). The non-detection for L5 does 
not mean we find that the bias is linear on all scales, only 
above ~ 4/i~^Mpc. Quantitatively, we find a Ax^ for the 
best-fit model with 62 free versus that with 62 = (lin- 
ear bias only) of (0.2, 32, 4.5) for Main-L5, LRG, and LRG- 
highz, respectively. That number for Main-L5 confirms the 
non-detection of non-linear bias, but the Ax^ for LRG sug- 
gests a > 3(T detection of non-linear bias, with a marginal de- 
tection for LRG-highz. As expected, the best-fit model with 
62 = has higher erg = 0.85, 0.83 for LRG and LRG-highz, 
to accommodate the increased clustering below lOft^^Mpc 
due to non-linear bias. 

As described in Sec. 12.3.41 we have included an arbi- 
trary galaxy-galaxy lensing signal calibration uncertainty 
with a Gaussian standard deviation of (0.04, 0.04, 0.05) for 
the three samples, assuming the calibration is perfectly cor- 
related between them. Given this prior, we find that the 
best fit is obtained with a calibration that is O.Ict from the 
expected value (this decrease in the signal is illustrated in 
Fig. [To] via an increase of the theoretical model). If we in- 
stead fix the lensing calibration without allowing any free- 
dom, the best-fit as changes by an amount that is well below 
our quoted error, and the errors become smaller by 20 per 
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Figure 10. Top left: Observed -RTg^ with best-fit signals from fits 1-3 in Table [S] for these fits, the data were fit separately for each 
sample. Data have been rebinned for easier viewing. Bottom left: Ratio of observed to best-fit signal from the top panel, using the original 
narrower binning as for the actual fit. Right: Same as left, but for clustering iJTgg. 
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Table 3. Fits for cosmological parameters described in Sec. 17.11 In each case, parameters that are fixed to a single value are in bold; 
those that arc fit are shown with 68 per cent confidence limits after marginalising over all other fitted parameters. 



cent. These findings suggest that systematic uncertainty due 
to uncertain lensing calibration is not completely negligible, 
but does not dominate our error budget. 

We have also carried out these fits with the off-diagonal 
elements of the covariance matrix set to zero for both lensing 
and clustering for all samples, in order to test how sensitive 
we are to the treatment of off-diagonal elements. We find 
that the best-fitting erg changes by 0.01 (0.2cr), and the size 
of the error regions becomes smaller by 15 per cent. Thus, 
our results are relatively stable to inclusion of correlations 
between radial bins. 

In our phenomenological bias model described in 
Sec. 12.21 we have included a relation 62, gm = &2,gg — A and 
fixed A = 0.2. To test how sensitive we are to this choice 
(which better fits the data in simulations than A = 0), we 
have redone these fits with A = 0. Typically, we find that the 
resulting as shifts by +0.027, or half of our combined statis- 
tical and systematic error. However, we are slightly sensitive 
to this change: with the same number of degrees of freedom, 
the fit with A = has a higher x^, with Ax^ = 3. This 
finding justifies our choice of A = 0.2 for all cosmological 
fits. 

Next, we have allowed the matter density Q.m to vary 



(fit 5 in Table O using all three samples together) . In this 
case, there is a classic degeneracy between ag and Q.m in the 
lensing data, with higher ag requiring a lower Q.m to fit the 
data. Because of this degeneracy, the allowed range of (Tg be- 
comes broader, erg = 0.76±0.08 when marginalised over fi^, 
the bias parameters, and the lensing calibration (still a 10 
per cent constraint on erg even with this degeneracy). When 
marginalising over erg, we can constrain 0.^ = 0.2mto°oll- 
The resulting 2D contour plot for these two parameters is 
shown in Fig. 1131 The best-constrained parameter combi- 
nation, which is a better illustration of our overall S/N , is 
as,{Q.m/Q.2^)°'^'^ = 0.80±0.05 (Icr, stat. + sys.), representing 
a 6 per cent uncertainty in the amplitude of matter fluctu- 
ations. We do not show the best-fltting signal for this case, 
because it differs very little from that shown in Fig. 1101 



7.2 Combination with other data 

We combine these data with WMAP7 CMB results 
(|Komatsu et al.ll2011^ ■ anticipating significant benefit from 
a tight prior on the amplitude of matter fluctuations 
at early times. For the flrst combined fit, we vary 
(A, Wdm, Wb, 6', T, ris) (see definitions at start of Sec. 17. ip . 
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Figure 11. Best- fitting galaxy bias parameters for Main-L5 (top), LRG (middle), and LRG-highz (bottom). Left: The posterior proba- 
bility distribution for the large-scale galaxy bias b, marginalised over other parameters including the non-linear bias parameter 62. The 
solid black lines show the results for fits 1-3, fitting for each sample separately; the dashed red lines are the results from fit 4, jointly 
fitting all samples. Right: Contour plots for the large-scale bias b versus non-linear bias parameter 62 (1, 2, and So"). The black line 
contours show the results from fits 1-3, fitting for each sample separately; the coloured contours show the results from fit 4, jointly fitting 
all samples. 



the six bias parameters for the three galaxy samples, and 
the lensing calibration with its (4, 5, 5) per cent Gaussian 
prior - with fixed w^c = — 1 and Qk ~ 0. We include lens- 
ing of the CMB and m arginalise over an SZ template as in 
iKomatsu et al.l (120111 ). All results in this section that use 
our data include our (t| prior from Appendix |X1 but that 
prior was not included for the analyses that use WMAP7 
data alone. 

Fig. [14] shows 2D parameter contours with WMAP7 
alone and combined with these data (the Sim vs. erg con- 
tour was already shown in Fig. I13[) . We can see that since 
the degeneracy direction for as vs. for the CMB data 
is orthogonal to that for lensing data (Fig. I13p . adding our 
data to the WMAP7 data roughly halves the size of the al- 
lowed regions in parameter space. However, the top panel 
of Fig. [14] shows that the constraints on Us are, unsurpris- 
ingly, set exclusively by the WMAP7 data. The bottom 
panel shows that we also provide additional constraining 
power on Ho, through the stronger constraint on Qrn (and 



the fact that CMB constrains a different parameter com- 
bination, Qrnh^)- The resulting ID probability distributions 
for (hs, Us, ^m. Ho), marginalised over other parameters, are 
shown in Fig. 1151 

Table [4] gives best-fitting parameters and their 68 per 
cent confidence intervals from these fits, for WMAP7 alone 
and for WMAP7 plus these data (fits 6 and 7). As shown 
in Figs. [14] and 1151 the parameters for which constraints 
improve significantly by combination of these datasets are 
as, fim, and Ho- In these combined fits, there is little 
tension between the datasets, and the best-fitting galaxy 
bias and lensing calibration bias parameters are largely un- 
changed from their values when fitting only to the SDSS 
data. As in the previous section, we can identify the best- 
fitting combination of erg and Q,m, which has changed 
from crg(nm/0.25)"-^^ = 0.80 ± 0.05 (our data alone) to 
crg(!:!m/0.25)"° " = 0.79 ± 0.02 (WMAP7 + our data) 

Next, we allow the equation of state of dark energy 
(wdc) to vary from a cosmological constant. Figs. 1161 and ll7l 
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Fit Data eg Q,m «s ^^o "^dc ^fc X] "^i^ 

km Mpc~l cV 



6 WMAP7 0.810 ±0.029 0.270Io;o27 0.965 ± 0.014 70.4 ±2.5 -1 

7 WMAP7+our data 0.796 ± 0.019 0.261 ±0.014 0.966 ± 0.013 71.1 ±1.5 -1 

8 WMAP7 0.83t° :['|' 0.26j;o J° 0.969 ± 0.014 72 ± 11 -1.05^0.30 

9 WMAP7+our data 0.82 ± 0.08 0.25lo°3 0.968 ± 0.014 73 ± 6 -1.07 ±0.20 

10 WMAP7 0-77to;o? 0.42j;j;;^J, 0.964 ± 0.015 57^^^ -0.9Atofi -0-03t^;^^ 

11 WMAP7+our data 0.81^0.08 ^■'^^-o'^ot 0.969 ± 0.014 72^^ -l-05to 24 0.00 ± 0.01 

12 WMAP7 0. 72^0 08 0.32j;J5;°^ 0.961 ± 0.016 66 ± 4 -1 < 1.1 

13 WMAP7+our data 0.76^0.05 0.29,^1'°^^ 0.968 ± 0.013 69 ± 2 -1 < 0.56 



Table 4. Fits for cosmological parameters using WMAP7 and our data as described in Sec. 17.21 In each case, parameters that are fixed 
to a single value are in bold; those that are fit are shown with 68 per cent confidence limits after marginalising over all other fitted 
parameters. The exception to this convention is the sum of neutrino masses, for which the 95 per cent upper limit is shown. 




Figure 12. Posterior distribution for erg marginalised over all 
other fit parameters, for fits 1-4 (Om fixed). 



show 2D contours and ID parameter distributions, respec- 
tively, and Table |4] gives best-fitting parameter constraints, 
again with WMAP7 alone (fit 8) and with our data included 
(fit 9). As in the ACDM case, our data do not provide sig- 
nificant additional constraining power on ris, but it does 
improve the constraints on all the other parameters, most 
dramatically on Wdo (because adding our data provides a 
constraint on the amplitude of matter fiuctuations at times 
well after dark energy has become important). 

Next, we relax the assumption of flatness (while still 
allowing the equation of state of dark energy to vary from 
a cosmological constant). Figs. [181 and 1191 show 2D and ID 
constraints, respectively, and Table [4] gives best-fitting pa- 
rameter constraints, again with WMAP7 alone (fit 10) and 
with our data included (fit 11). We can see that with these 
relaxed assumptions about Wde and Q.k, the posterior prob- 
ability distributions for O^, Hq, and Wde are very broad; 
our data play a crucial role in reducing the allowed region 
of parameter space for all parameters except . 



Figure 13. Open black lines show the contours in the trg vs. 
Om plane (fit 5) for our dataset, marginalising over all linear and 
non-linear bias parameters and lensing calibration. The contours 
that are shown are 1, 2, and 3(T. The nearly orthogonal open 
red contours are for WMAP7 (also fitting for Us, Ho, and other 
parameters in fiat ACDM as in Sec. I7.2( l. and the filled contours 
are for WMAP7 combined with our data. 



Finally, we revert to the assumption of flat ACDM, but 
we allow for the presence of massive neutrinos using the for- 
malism in Sec. 12.31 Fig. [20] shows ID posterior probability 
distributions and Table [4] gives best-fitting parameter con- 
straints, again with WMAP7 alone (fit 12) and with our data 
included (fit 13). In addition to tightening constraints on 
erg, Qrn, and Ho as before, our data help to rule out models 
with neutrino masses on the higher end of those allowed by 
WMAP7, such that the one-sided 95 per cent upper limit on 
neutrino mass goes down from 1.1 eV with WMAP7 alone to 
0.56 eV with WMAP7 and our data. In general, the presence 
of massive neutrinos can significantly broaden the allowed 
parameter space for Qrn and Ho from CMB alone, and our 
data rule out some of the allowed high-Qm or low-//o values. 
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Figure 15. ID probability distributions for fits using the 
WMAP7 data along with our new results for galaxy-galaxy lens- 
ing and galaxy clustering in SDSS, assuming flat ACDM as in 
Fig. 1141 In all cases, we have marginalised over the nuisance pa- 
rameters (bias and calibration) for our data, and over any cosmo- 
logical parameters not shown on the plot. The solid black lines 
are for WMAP7 alone; red dashed lines are for WMAP7 and our 
data together. 
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Figure 14. Contour plots of 2D probability distributions for fits 
using the WMAP7 data along with our new results for galaxy- 
galaxy lensing and galaxy clustering in SDSS, assuming flat 
ACDM. In all cases, we have marginalised over the nuisance pa- 
rameters (bias and calibration) for our data, and over any cosmo- 
logical parameters not shown on the plot. The black lines show 
1, 2, and 3a contours for WMAP7 alone. The contours shown 
in colour are WMAP7 and our data together. Top: erg vs. Us; 
bottom: erg vs. Hq. 



7.3 Comparison with previous work 

The results of the previous section illustrate that our data 
provide cosmological parameter constraints that are consis- 
tent with and complementary to those from WMAP7. Here, 
we compare our constraints with those from other cosmo- 
logical probes. 

We do not compare our results against those from a pure 
clustering analysis of the shape and amplitude of the galaxy 
power spectrum, due to systematic uncertainties in treat- 
ments of non-linear galaxy bias, redshift-space distortions, 
and other issues. However, it is valuable to compare against 
measurements of baryonic acoustic oscillations (BAO), a 
measure of the expansion history of the universe rather than 
the growth of structure, since it is significantly less prone to 
such uncertainties. The most recent measurement of BAO 
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Figure 16. 2D contour plot for flat uiCDM flts with the equation 
of state of dark energy allowed to vary from our fiducial value of 
— 1 (but assumed to be constant in time). Line and contour styles 
are as in Fig. 1141 
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Figure 17. ID probability distributions for fiat wCDM fits witfi 
the equation of state of dark energy allowed to vary from our 
fiducial value of —1 (but assumed to be constant in time), as in 
Fig. If 61 Line styles are as in Fig. 1151 



Figure 19. ID probability distributions for wCDM fits writhout 
the assumption of flatness, as in Fig. 1181 Line styles are as in 
Fig. [15] 



in the SPSS DR7 jMehta et al.l I2OI2I : IPadmanabhan et all 
|2012| : IXu et aL| [ 2012l ) represents the first us e of the 'recon- 
struction' technique (jEisenstein et al.] |2007| ) to reduce the 
effects of non-hnear evolution of the den sity field in smooth- 
ing the BAO peak. iMehta et al.l (|2012l ) demonstrates that 
when combining the BAO peak position with CMB data, 
they find fi^ = 0.280 ±0.014 (in the context of flat ACDM). 
This result is fully consistent with our data when fitting for 
as and Q.m, which yielded fi™ = 0.257+°;o34 (Sec. [73]). This 
consistency is non-trivial, given that the BAO results use 
an identical sample to make a fully geometric constraint on 
cosmology, whereas we measure the amplitude of clustering 
well below the BAO peak. 

We can also compare against the BAO results from the 
B aryon Oscillation Spec troscopic Survey (BOSS) presented 
in lAnderson et al] l|2012l '). Assuming flat ACDM, when using 
WMAP7 and the BAO results in that paper alone, they find 
Q.rr, = 0.298±0.017 and Ho = 71.0±2.6km s"^ Mpc'^ The 
Q,m constraint is ~ Icr above ours, consistent within the 
errors. 

We also compare our results against those from other 
lensing analyses, particularly cosmic shear. First, we com- 
pare against those fro m the COSMOS sur vey, including the 
or iginal analysis fromlMassev et all (|2007l ) an d a re-analysis 
in ISchrabback et al.l (|2010l ). The results from iMassev et al] 
l|2007^ used a 3D analysis to infer 0-8(^^/0.3)°'"' = 
0.866lo:o68 (68 per cent CL, stat. + sys.). We can com- 
pare this result against our result when fitting for irg and 
fim, cr8(J7™/0.25)° " = 0.80 ±0.05 marginalised over nui- 
sance parameters. The COSMOS results are ~ l.Qa above 
ours, giving a higher am plitude of clustering. The 3D COS- 
MOS lensing analysis in ISchrabback et al.l (|2010l ) gives, for 
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flat ACDM, a value cr8(nm/0.3)° 



: 0.75±0.08, consistent 



Figure 20. ID probability distributions for flat ACDM fits with 
massive neutrinos. Line styles are as in Fig. 1151 



with our results at the ~ 0.2o- level. Part of the reason for 
the lower quoted clustering amplitude in Schrabbac k et al] 
(2010) is a different treatment of the non-linear power spec- 
trum (m ore consistent with o urs) : if they use the same treat- 
ment as iMassev et al.l l|2007l ). they find 0.79 ± 0.09, higher 
by 5 per cent. Other differences in clustering amplitude be- 
tween the two COSMOS results could come from the dif- 
ferent treatment of PSF estimation, charge-transfer ineffi- 
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Figure 18. 2D contour plot for liiCDM fits without the assumption of flatness. Line and contour styles are as in Fig. 1141 



ciency, the availability of more photometric data to improve 
the photometric redshifts, or diflterences in analysis meth- 
ods (scales used and so on). In short, the COSMOS lensing 
results are consistent with ours, with the exact comparison 
depending on the method of analysis and the treatment of 
systematic errors. 

We can also compare against cosmic shear results from 
stripe 82 of the SDSS itself. There are two such results 
that use largely ind e pende nt a nalysis methods o n the same 



ttiat use largely mcl e pende nt a nalysis methods o n the same 
d ata, bvlHuff et al.l (120111 ) and iLin et al] (|201ll ). The work 
in lHuff et al. (2011) used the same PSF correction technique 
as we have used, and also the sam e simulation method rely - 
ing on space-based training data l|Mandelbaum et al.ll201 
to calibrate the shape measurements, so its systematic er- 



rors may not be fully independent from ours. However, given 
that the area of stripe 82 is ~ 3 per cent of the area used 
here, we can consider those results to be statisti cally inde- 
ende nt of ours. With fixed D,m close to our value. iHuff et al.l 
20111) find a relatively low amplitude of matter fluctua- 
tions, erg = 0.636^0 {54, which can be compared with our 
0.80 ± 0.05. Assuming completely independent errors, this 
represents a 1.4a discrepancy, which is not statistically sig- 
nificant mainly due to the small size of stripe 82 and the 
limited sour ce number densit y due to the SDSS seeing. Com- 
paring with lLin et"aLl l|201ll) . they find for flat ACDM that 
asi^m^ = 0.252lg;[Jg. For our value of fi™ = 0.25, that con- 
straint becomes erg = 0.665101371 quite similar to that from 
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iHuff et al.1 (|201ll ). In both cases we therefore see a shght 
tension with our results, but only at the ~ lAa level. 

We also compare against the cosmic shear results 
from the Canada- France-Hawaii Tel escope Legacy Surve y 
(CFHTLS) from several years ago (|Hoekstra et al.ll2006h . 
Using 22 deg^ for fixed Um = 0.3, they find as = 0.85±0.06. 
Rescaling this to fim = 0.25 using us oc STi-" '', their re- 
suh is 0.95 ± 0.07. Re sults for the CFHTLS Deep survey 
(|Semboloni et al.1 120061 ') were similar. While this value is 
1.8(7 above our result at ^Im = 0.25, erg = 0.80 ± 0.05, an- 
other paper b ased on an incre ased area working only in the 
linear regime (|Fu et al.l liooi ) reported asinm/0.25)°-''''^ = 
0.837 ± 0.084, more consistent with our results. In short, it 
is unclear whether there is real tension here; a comparison 
with the full CFHTLS dataset will be more informative. 



8 DISCUSSION 

We have used updated measurements of galaxy-galaxy weak 
lensing and galax;y clustering for several samples of spec- 
troscopic galaxies in the SDSS DR7 to place competitive 
constraints on the amplitude of matter fiuctuations and, 
by combining with WMAP7 data, the growth of structure 
with time. The novelty in comparison with previous lensing 
cosmology constraints is that we have used galaxy-galaxy 
lensing (a cross-correlation of lens galaxy positions and the 
background shear field) rather than cosmic shear (the auto- 
correlation of the shear field). From a statistical perspective, 
the former is more detectable in low-redshift surveys such as 
SDSS; however, even at higher redshift, the galaxy-galaxy 
lensing typically has a lower systematic error budget, be- 
cause the use of cross-correlations allows us to more easily 
remove several systematic errors (intrinsic alignments, addi- 
tive shear systematics) that plague cosmic shear. To avoid 
contamination from the smallest scales, where there are un- 
certainties in the galaxy-mass cross-correlation due to the 
way that galaxies populate halos, we have used the annu- 
lar differential surface density (ADSD) statistic T, which 
strictly removes contributions below some scale -Ro (cho- 
sen based on comparison with simulations to be 4/i~^Mpc 
for the clustering analysis, and 2/i~^Mpc for the lensing 
analysis). We apply this approach to three different non- 
overlapping samples extracted from SDSS DR7: an interme- 
diate redshift (0.16 < z < 0.36) LRG sample, high redshift 
(0.36 < z < 0.47) LRG sample, and a low redshift sample 
{z < 0.155) with a fainter absolute magnitude limit and no 
colour selection. We have opted to model the signals using 
the non-linear matter power spectrum along with a pertur- 
bation theory-based model for the non-linear galaxy bias, 
containing parameters over which we marginalise. We see 
clear evidence for a scale-dependent bias in our LRG sam- 
ple with large-scale bias around 2, while there is no evidence 
for the scale-dependent bias for the lower luminosity sample 
with large-scale bias around 1.25. This trend is consistent 
with theoretical expectations based on simulations and an- 
alytic predictions (Balda,uf et al. 2012). 

Using our data and fixing all cosmological parameters 
except for as and Qm, we find cr8(n„/0.25)" " = 0.80±0.05 
{la, stat. -I- sys.) after marginalising over the galaxy bias 
parameters and a nuisance parameter for the lensing cali- 
bration. This result is highly consistent with that from the 



WMAP7 CMB analysis, and with many other cosmologi- 
cal measurements as discussed in Sec. 17.31 The 6 per cent 
errorbar, including both statistical and systematic contribu- 
tions, indicates that the SDSS is a quite powerful survey for 
weak lensing cosmology at low redshift (in the context of 
other extant lensing datasets). Moreover, given its low effec- 
tive redshift of ~ 0.27, we imagine future benefits from the 
combination with other lensing measurements that typically 
have higher effective redshifts. 

When we include WMAP7 data in the analysis, we find 
that for flat ACDM, we are able to provide significant addi- 
tional constraining power on erg, Q,rn, and Hq due to orthog- 
onal degeneracy directions, effectively halving the allowed 
region in parameter space; we do not provide significant ad- 
ditional constraining power on Us . When we allow the equa- 
tion of state of dark energy Wdc to vary from —1 (while still 
assuming it is constant in time), when we further allow the 
possibility of curvature, or when we include massive neutri- 
nos in the context of fiat ACDM, we likewise find that our 
low-redshift constraint on the amplitude of matter fluctua- 
tions is crucial for reducing major parameter degeneracies. 
It will be interesting to combine these results with a low- 
redshift constraint on the expansion history of the universe, 
such as from BAO; we defer this exercise to future work. 

We emphasise that these results represent an entirely 
new opportunity for the field of lensing to constrain cosmo- 
logical parameters in a way that is largely independent of the 
details of how galaxies populate dark matter halos, but also 
less sensitive than cosmic shear to several important obser- 
vational and astrophysical systematic errors. Among these 
are all additive contributions such as telescope or atmo- 
sphere effects that induce shear-shear correlations but not 
shear-galaxy correlations. Intrinsic alignments also do not 
contribute to shear-galaxy correlations as long as the red- 
shift separation between lenses and sources, using photomet- 
ric redshift information, is effective. In contrast, the domi- 
nant intrinsic alignment contribution to shear-shear corre- 
lations, induced by correlations between sheared galaxies in 
the background and intrinsically-aligned galaxies in the fore- 
ground, cannot be el iminated simply by us ing photometric 
redshift information l|Hirata fc Seliak|[2004l '). 

In addition to the intrinsic value of these cosmo- 
logical constraints in and of themselves, this work is 
a proof of concept for this analysis technique for the 
next generation of large, wide-field surveys that will 
carry out lensing meas uremen t s, suc h as Hyper Suprime- 
Cam (HSd3, Mivaza ki et alj I2006D. Dark Energy Sur - 
vey (DE^^T frhe Dark Energy Survey Collaboratiorji2005l '). 
the Kilo-Degree Survey (kidC!|), the Panoramic Survey 
Telescope and Ra pid Response System ( Pan- S TARRa ^ T 
iKaiser et al.ll2O10l '): and even more ambitious programmes 
such as the Large Synoptic Survey T elescope (LSSTpl. 
ILSST Science CoUaborations et aLll2009l '). EuclicQ, and the 



^® http://ww.naoj . org/Projects/HSC/index .html 
https : //www . darkenergysurvey . org/ 
http : / /www . astro-wise . org/pro j ects/KIDS/ 
http : //pan-Starrs . if a.hawaii . edu/public/ 
http://www.lsst.org/lsst 

http: //sci . esa. int/science-e/www/area/ index . cf m?f areaid 
=102 
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Wide-Field Infrared Survey Telescope (WFIRSlH). The 
ability to make cosmological measurements with galaxy- 
galaxy lensing rather than cosmic shear is particularly im- 
portant for making use of early data from upcoming surveys, 
when additive shear systematics will be less well-understood. 
We expect that this approach will yield results that are com- 
petitive and complementary to shear-shear analysis for these 
future surveys as well. 
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APPENDIX A: CALIBRATION OF THE 
METHOD 

As a basic test, we use the data from the simulated LRG 
sample (Sec. 12. ip and check whether we can accurately re- 
cover the input cosmology using the analysis framework 
from Sec. 12.31 Such tests with various i?o will allow us to 
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determine the safest choice of Rq to minimise systematic 
error without excessively increasing the statistical errors. 

For this test, we fixed all cosmological parameters be- 
sides as , and varied _Ro , making sure that we can recover the 
true as = 0.8 for the simulationf^. While doing this test, we 
allowed the lensing calibration and the bias parameters to 
be free parameters. Given that moving Ro above 2/i~^Mpc 
for the galaxy-galaxy lensing severely impacts the S/N, we 
consider only this value of Ro for the lensing, but vary Ro for 
the clustering from 2 to 6/i~^Mpc, in steps of 2/i~^Mpc. This 
procedure is also theoretically motivated since, as pointed 
out in Sec. 12.41 the clustering signal should be more strongly 
sensitive than the g-g lensing signal to the details of how 
galaxies occupy dark matter halos on these scales. Shape 
noise was not added to the simulated LRG data, but we 
used the covariance matrix from the data, so biases that 
are present at the level of Acts ^ 0.01 are true biases rather 
than noise fluctuations. For -Ro.gg = 2, 4, 6/i~^Mpc, the best- 
fitting CT8 = 0.763,0.795,0.792 (these are the global best- 
fit values, not marginalised over nuisance parameters). The 
first value, for -Ro,gg = 2/i~^Mpc, indicates a statistically- 
significant bias in as- However, the results for larger _Ro are 
consistent with no bias, so we adopt _Ro = 2/i~'^Mpc and 
4/i~^Mpc for the g-g lensing and galaxy clustering, respec- 
tively. 

We can also check the effect of priors and marginalisa- 
tion over lensing calibration bias and galaxy bias parame- 
ters. After marginalisation, the median of the likelihood is 
0.78, which difi'ers from the input value of 0.8. Since our 
best fit model (maximum likelihood) is at the position of 
the input model, the discrepancy between the median and 
the input value must be due to the effect of the prior. The 
median is the standard value quoted in the Markov Chain 
Monte Carlo (MCMC) analyses, since it can be robustly 
estimated. It is also invariant under a monotonic transfor- 
mation of the variable, i.e. the median of as is the same 
as the median of ct| assuming the same prior, which is be- 
lieved to be useful since there is no good reason a priori to 
use the linear fiuctuation amplitude cts as opposed to the 
correlation function amplitude ct|. However, this does not 
mean that the median is invariant under the change of the 
prior, i.e. if we start with a uniform prior on as we obtain a 
different median than if we use a uniform prior on ct|. The 
only number that is prior-independent is the maximum like- 
lihood value, which is unstable in multi-dimensional MCMC 
analyses, since the likelihood surface is shallow and MCMC 
has difficulties finding the absolute maximum with a finite 
number of steps. It is nevertheless convenient that the main 
reported number agrees with the input value in simulations. 
Since a prior that is uniform in cts does not result in this 
property, and since there is nothing especially natural about 
that choice of the prior, we will empirically choose a prior 
such that the median cts agrees with the input value in simu- 
lations, and then report the median determined in a similar 

We could have fit jointly for erg and Qm- However, since those 
two cosmological parameters are strongly degenerate, we could 
easily (due to very small noise fluctuations in the simulations) 
be driven anywhere along that degeneracy, without it being a 
meaningful deviation. We therefore check that with Qrn fixed we 
can infer the correct erg , under the assumption that this shows we 
can infer the proper degenerate combination of f!,n and erg. 



way in the data. If we apply a flat prior on CTg, the median 
increases to cts = 0.80, and since this agrees with the input 
value we adopt this prior for the rest of the paper. Note that 
as we add more data, such as additional SDSS datasets (we 
use three samples, not just LRCs as in the tests in this sec- 
tion) and WMAP, the effect of the prior is diminished and 
we converge on the true value. 

The plots of simulated T are shown in Fig.[l] and we see 
that the simulated data and the best model agree reasonably 
well both for galaxy-galaxy lensing and galaxy clustering, 
with b2 = 0.25 as the best fit value. Note that the same 
value of 62 also fits the cross-correlation coefficient for ^, as 
seen in the top panel of Fig. [1] 

We also note that in the real data, when carrying out 
the fitting procedure, we see a qualitatively similar trend 
in CTg(i?o,gg) increasing with Ro, however the change in cts 
is larger in amplitude (but consistent with that given here 
within the errors). As a final sanity check of our procedure 
above, we carry out the analysis using another HOD sample, 
constructed from simulations with cts — 0.90 and selected to 
mimic a higher-luminosity and higher-mass sample. For this 
HOD sample, with our adopted values of Ro for g-g lensing 
and galaxy clustering and our prior on cts, we find best- 
fitting and median cts = 0.90 and 0.92. The former number 
is completely consistent with the input cosmological model, 
reassuring us that our choice of Ro is robust to significant 
changes in the galaxy mass and input cosmology. The latter 
number suggests that for large values of cts, the adopted 
prior will cause a bias that is ~ 40 per cent of our statistical 
errors. However, we see no evidence for such a high cts, and 
so this worst-case scenario is not an issue in practise. 



APPENDIX B: WEAK LENSING 
SYSTEMATICS TESTS 

In this Appendix, we present several tests of systematic error 
in the lensing signals, focusing exclusively on those tests that 
were not done for the same shape catalogue by Rll. 

Bl Calculation of Tgm 

Since we use Tgm for cosmological parameter constraints, 
and it is a derived quantity that relies on determination of 
AE(i?o), here we present tests illustrating the accuracy of 
that determination. 

As stated in Sec. 14.11 we determine AE(_Ro) using 
power-law fits over a range of scales on which AE appears 
consistent with a power-law, including Ro itself to avoid ex- 
trapolation. In Fig. IBll we show (for all three samples) the 
observed AE(7?) divided by the best-fit power-law for the 
chosen range of scales which are indicated by vertical lines. 
This power-law is determined from a fit to the jackknife 
mean AE, weighted by the inverse variance (assumed to be 
diagonal, which is appropriate for these scales). Ideally, this 
ratio of observed signal to power-law would be consistent 
with one for all scales between the vertical lines. In addi- 
tion, we show the jackknife mean value of AE(7?o) and its 
jackknife errorbar, again with respect to the best-fit power- 
law for the chosen range of scales and therefore with an ideal 
value of 1. It is clear that the power- law approximation is 
valid on the range of scales used, but the observed signal 
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Figure Bl. Lensing signals for the three samples (as labelled 
on the plot) divided by the best-fit power-law using the range 
of scales indicated by vertical solid lines. The horizontal dashed 
line shows the ideal value AE/ASpow = 1. Also, a single point 
shown as an X with its own errorbar at Ro = 2/i^^Mpc shows the 
jackknife mean value that was used for AS(i?o)i and its jackknife 
error. 



deviates from it strongly outside of that range. This power- 
law fit is used only for empirical determination of AE(i?o), 
not for any other purpose. 

The best-fit power-laws and ranges of scales used are de- 
fined as AE = AEo(7?/i?o)°' for R 

pow,min ow,max , 

where AEo is in units of /iMq/pc^. Given this definition, the 
power- laws that went into Fig. lBll are described in Table iBl] 



B2 Correction for physically associated sources 

The boost factors (Eq. I33|l that implicitly went into those 
signals are shown in FiE. IB2l As shown, the correction is very 
small (~ per cent level) on the scales used for this measure- 
ment, > 2/i"^Mpc. The size of the errorbars indicates that 
there are ~ 1 per cent-level density fluctuations in the real 
catalogue on large scales, perhaps due to dependence of the 
lens number density on systematics that are not properly re- 
produced in the random catalogues. This does not bias the 
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Figure B2. Boost factors B(R) — 1 as a function of separation 
for the three samples used in this analysis, as labeled on the plot. 
The vertical solid line at i? = 2h~^Mpc indicates the physical 
scales used for the cosmological constraints in this paper. 



lensing signal, it simply acts as a minor contribution to the 
statistical error budget, subdominant to shape noise. 



B3 Intrinsic alignments 

We must consider the contribution of intrinsic alignments 
to our measurement. In principle, they should only con- 
tribute if photo-2; errors scatter galaxies that are actually 
physically-associated with our lenses into the source sample, 
and if those physically-associated galaxies have a tendency 
to align radially or tangentially with respect to our lenses. In 
principle, we expect a radial alignment for elliptical galax- 
ies, whi ch shou ld follow the linear alignment model to some 
extent llHirata fc Scliak 2004): this alignment has been ob- 
served (e.g.. iHirata et al.ll2007r ) but primarily for bright red 
galaxies, rather than the faint ones that we use as sources in 
our lensing analysis. Disk galaxies are expected to align in a 
way that relates their disk angular momentum to the dark 
matter halo angular momentum, which leads to essentially 
no net r adial or tangential ali gnment with respect to lens 
galaxies l|Hirata fc Seliakll2004l '). 

Thus, our approach here is to assume that blue sources, 
which contribute ~ 70 per cent of the weight to the lensing 
measurement, contribute zero intrinsic alignment, and red 
sources contribute some amount that we must empirically 
constrain. A method for carrying o ut th is empirical con- 
straint is presented in lBlazek et al.l l|2012l ): it relies on cal- 
culating shears of sources selected in different photo-z bins 
with respect to the lenses. That work shows that using this 
method, we can constrain intrinsic alignment contamination 
of the galaxy-galaxy lensing signal at 10/i~^Mpc (the scale 
which roughly dominates our lensing constraints) to < 2 per 
cent, for our LRG lens sample (which has the lowest statis- 
tical error and therefore is the most demanding in terms 
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of systematics) . Assuming that ~ 30 per cent of the source 
sample has at most a 2 per cent intrinsic ahgnment contam- 
ination, this translates to < 0.6 per cent intrinsic alignment 
contamination of our lensing signals. This is far subdomi- 
nant to our other assumed systematic errors (e.g., 4-5 per 
cent for lensing calibration) , so we assume it is negligible for 
the purpose of these cosmological constraints. 

B4 Ratio test 



As described in Sec. 14. 2. Tl we apply a number of corrections 
for known sources of calibration bias such as photo- 2; error. 
To support the claim that we understand these effects well, 
we carry out ratio tests of the signals computed using the 
same lens sample and different source samples. After cor- 
recting for known calibration biases that are different for 
each sample, these ratios should simply be consistent with 
one within the errors (which are typically ±0.06). 

There are several ratio tests that we can carry out, each 
of which is sensitive to our understanding of several effects. 
For example, when dividing the sample into red and blue 
galaxies, we are most sensitive to differences in photo-2 er- 
rors and to differences in the intrinsic ellipticity distribution 
of the red and blue galaxy populations. When dividing into 
faint versus bright galaxies, we are most sensitive to noise 
rectification bias and photo-z errors in the former. When 
dividing into well- vs. poorly-resolved galaxies, we are most 
sensitive to selection biases that couple the shear to the ap- 
parent size. When dividing into all sources vs. only those 
with Zphot > 0.45, we are most sensitive to photo-z errors in 
the former (the photo-z errors in the latter are less impor- 
tant because of the larger redshift separation between lenses 
and sources). Therefore, an ability to achieve a ratio of one 
in each of these four cases would suggest that we understand 
the different systematic errors that are affecting the lensing 
signals in each case well enough to constrain cosmology at 
the ~ 5 per cent level. 

Our findings, using the LRG lens sample, is that in three 
cases the ratio is within la of one. In the final case, the split 
into red versus blue sources, the ratio is 0.92 ± 0.06; thus 
the discrepancy is only a little more than la, and since we 
expect such a case to turn up after doing several ratio tests, 
we conclude that the ratio tests do not indicate any ma- 
jor misunderstanding of the predominant systematic errors 
affecting the lensing calibration. 

B5 Large-scale systematic shear 

In this section, we pre sent tests of the large-scale shear sys- 
tematics. As noted in iMandelbaum et al.l l|2005h . the pres- 
ence of a coherent PSF ellipticity along the scan direction 
in the SDSS data can result in a non-zero tangential shear 
on large scales, where lens-source pairs may be lost due to 
survey edge effects. There, we corrected for this effect by 
subtracting off the g-g lensing signal around random points, 
but also noted that if the systematic shear correlates in some 
way with the lens number density, this correction may not 
be sufficient. Thus, we must test the accuracy of this pro- 
cedure when using large-scale lensing signals to constrain 
cosmology. 

In Fig. IB3I we show the g-g lensing signal around ran- 
dom points for our three lens samples. As shown, it becomes 
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Figure B3. The black points and lines show —RAH^^^^, a mea- 
sure of systematic shear in the source catalogue. Each panel is 
a different sample (Main L5, LRG, LRG-highz from top to bot- 
tom). For reference, the actual signal that we use for science is 
also shown as the red dashed lines. 



significantly inconsistent with zero for scales above 10, 20, 
and 15/i~^Mpc for Main-L5, LRG, and high-z LRG, respec- 
tively. The reason for the different scales is that it is a sys- 
tematic associated with an angular scale, which becomes dif- 
ferent physical scales when we convert angular separation 6 
to R at the different typical lens redshifts. However, its mag- 
nitude also depends on the apparent magnitude and resolu- 
tion factors of sources used, as we have confirmed explicitly 
by dividing our source sample by source propertiej^. These 
two effects go in competing directions, driving the typical 
scale of the systematic error to increasing R as z increases, 
then at a certain point above z ~ 0.3, the typical scale starts 
to decrease again. For reference, this figure also shows the 
observed signal around real lenses in the same form (TiAS), 
to show at which point the systematic correction is com- 
parable in size. Generally, on the maximum scale used for 
science, the systematic correction ranges from 1/3 to 1/2 of 
the real signal, and is 1-2 times larger than the statistical 
error. 

Given the significance of this systematics signal on the 
largest scales used for our measurements, we must assess 
whether the assumptions behind the correction for it are cor- 
rect (c.f. Sec. l4.2|) . In order to do so, we rely on the following 
facts: First, if there is some coherent systematic shear, then 
depending on the survey geometry, it can also show up as 
a signal in the other shear component 7x on large scales. 
Second, there is no gravitational lensing signal in that shear 
component. Thus, we measure AEx around real lenses and 



This result is expected, because systematic errors in determi 
nation of the shear generically depend on both the S/N and ap- 



arent size of the source compared to the PSF; e.g.. iBridle et ahl 
2OIOI') . 
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Figure B4. Top; The black solid and red dashed lines show 
-RASx and iJAEx randi respectively, for the Main L5 sample. 
Bottom: For the scales that have a significant systematic shear 
signal, we show the ratio AEx,rand/^Sx • 
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Table CI. Best-fitting power-law functions (as defined in 
Sec. IClll to Wgg{R) for a limited range of scales, used to esti- 
mate uigg(iJo) and therefore Tgg. 



suiting systematic shear signal was consistent with zero, and 
the value of AE was unchanged on average on small scales, 
and decreased at the expected level on the largest scales. 
However, we merely present this test to validate our cor- 
rection scheme; we do not use the signal with survey edges 
removed for science, because the statistical errors on a few 
/i^^Mpc scales increase significantly (20 per cent). 

Finally, we note that a non- negligible fraction of this 
systematic shear signal (~ 30-40 per cent) is due to an error 
in the SDSS PSF model identified in Rll, which resulted in 
all PSF models in the r band in one of the camera columns 
having some spurious ellipticity. The remainder is due to 
the inadequacy of the adopted PSF correction method at 
removing the PSF ellipticity from the galaxy shapes. 



AEx.rand arouud random points, and check to see whether 
they are consistent. An inconsistency would call into ques- 
tion the validity of this correction using random points. 

As shown in Fig. IB4I for one of the lens samples, we 
find that there is a significant signal in the x shear com- 
ponent on the same scales as for the + component, and 
AEx.rand ~ 1.25AEx (actually it is ~ 1.2AEx on this Fig- 
ure; 1.25 is the result if we average over all the lens samples 
to reduce the noise). Taken at face value, this suggests that 
the systematic shear correlates with some factor that deter- 
mines the effective lens number density, and that correlation 
is not properly taken into account by our procedure for pro- 
ducing random catalogues. This finding is also consistent 
with the results in Sec. IB2I There is likely a simple physi- 
cal explanation for this issue related to the data processing; 
for example, very slight correlations of lens selection proba- 
bility with the PSF FWHM, extinction, or sky level should 
also correlate with fluctuations in the systematic shear in 
the sources. 

We consider the implications for the LRG sample, as an 
example. That sample has AE ~ 0.30 ± 0.03 and AErand ~ 
—0.10 at 50^~^Mpc. If we assume that AE^and is overesti- 
mating the true systematic shear by a factor of 1.25, then 
that means we should have been subtracting —0.08, not 
—0.10, from the original signal to get our final AE. Thus, 
our original AE should have been 0.28 ± 0.03, a shift of 
(2/3)cr. The size of this correction compared to the statisti- 
cal error suggests that we should impose the correction by 
simply dividing AE^and by 1.25 before subtracting it to get 
our estimate of the lensing signal. Indeed, we applied cor- 
rections in this way to all signals used for science and shown 
in all plots in this paper. 

As an additional test, we computed the signal while ex- 
cluding all LRGs within 60/i~^Mpc of a survey edge. The re- 



APPENDIX C: GALAXY CLUSTERING 
SYSTEMATICS TESTS 

CI Calculation of Tgg 

Since we use Tgg for our cosmological parameter constraints, 
and it is a derived quantity that relies on determination of 
Wgg{Ro), here we present tests illustrating the accuracy of 
that determination. 

As for the lensing signal AE, we determine Wgg(_Ro) 
using power-law fits over a range of scales including Ro on 
which Wgg appears consistent with a power-law. In Fig. ICll 
we show (for all three samples) the observed Wgg{R) divided 
by the best-fit power-law for the chosen range of scales which 
are indicated by vertical lines. This power-law is determined 
from a fit to the jackknife mean Wgg, weighted by the inverse 
variance (assumed to be diagonal on these scales). Ideally, 
this ratio of observed signal to best-fit power law would be 
consistent with one between the vertical lines. It is clear 
that the power-law approximation is valid on the range of 
scales used, but the observed signal deviates from it strongly 
outside of that range, which is what we expect for ACDM 
and typical models of scale-dependent bias. 

The best-fit power-laws and ranges of scales used are 
defined as uigg = wo {R/Ro f for Rp 

ow,min < R< Rp ow.max , 
where wo is in units of /i~^Mpc. Given this definition, the 
power-laws that went into Fig. ICll are given in Table ICll 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy-galaxy lensing cosmology 33 



2 
1.5 

1 

0.5 



2 

|l.5 

1 


2 

1.5 

1 

0.5 




i — I I I I I 



n 



'i i 



Main— L5 



! i 



i i 



H — I I I I I 



ii^ t , • I 1 i • • , 



H — I — li ne 



LRG 
" « 1 



» i 



i i 



H — I I I I I 



H — I — li ne 



LRG-highz 



10 



R [h-iMpc] 

Figure CI. Clustering signals for the three samples (as labelled 
on the plot) divided by the best-fit power-law using the range of 
scales indicated by vertical solid lines. The horizontal dashed line 
shows the ideal value ifgg/wgg.pow = 1- 
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